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RESEARCH  GOALS 

The  goal  of  this  project  was  to  develop  a  spiral-wave  front  beacon  for  UUV  navigation. 
Over  the  course  of  this  project,  the  spiral  beacon  was  brought  from  drawing  board  to 
prototype  in  collaboration  with  Ben  Dzikowicz  at  NRL.  The  performance  of  the  beacon 
has  been  evaluated  and  tested  both  in  a  series  of  field  experiments  and  through  numerical 
simulations.  The  final  stage  of  this  project  brought  the  spiral  beacon  from  a  technical 
readiness  level  (TRL)  of  4,  component  and/or  breadboard  validation  in  laboratory 
environment,  to  TRL  6,  system/subsystem  model  or  prototype  demonstration  in  a 
relevant  environment. 

FY08-FY09 


The  beacon  required  the  development  of  a  novel  transducer  capable  of  generating  a  spiral 
wave  front.  The  goal  of  the  first  two  years  of  this  project  was  to  develop  and  test  simple 
prototypes  of  the  spiral  wave  front  transducer  in  the  laboratory.  These  prototypes  were  to 
be  evaluated  both  numerically  and  analytically. 

Accomplishments 

Two  spiral  wave  front  beacon  prototypes  were  developed  and  tested  in  a  water  tank  at  the 
Naval  Surface  Warfare  Center-Panama  City  Division  (NSWC-PCD,  Panama  City,  FL, 
USA)  in  the  summer  of  FY09.  These  prototypes  were  also  evaluated  using  simple 
analytic  models  and  numerically  using  finite  element  methods.  The  results  of  the 
modeling  efforts  were  published  in  a  paper  titled,  “A  spiral  wave  front  beacon  for 
underwater  navigation:  Basic  concept  and  modeling  [1].” 


FY10-FY11 


The  goal  in  the  second  two  years  of  funding  was  to  develop  more  robust  versions  of  the 
spiral  wave  front  beacons  and  evaluate  the  performance  in  a  series  of  field  experiments. 
These  field  tests  were  also  accompanied  by  the  development  of  modeling  techniques  to 
assess  the  performance  of  the  beacon  in  operational  environments. 

Accomplishments 

In  FY10,  Dr.  Thomas  Howarth  and  Kim  Benjamin  at  NAVSEA  Newport  Division  in 
Newport,  RI  were  contracted  to  produce  two  new  spiral  wave  front  beacon  prototypes. 
The  results  of  laboratory  evaluations  of  these  prototypes  were  published  in  a  paper  titled, 
“A  spiral  wave  front  beacon  for  underwater  navigation:  Transducer  prototypes  and 
testing  [2].” 

These  prototypes  were  deployed  in  field  tests  at  the  NSWC  PCD  Acoustic  Test  Facility  in 
FY10  and  at  the  Navy’s  Dodge  Pond  facility  operated  by  the  Naval  Undersea  Warfare 
Center-Division  Newport  (NUWC-Newport,  Newport,  RI,  USA)  in  FY 11.  In  both  tests, 
the  beacon  was  secured  to  a  stable  platform  while  an  unmanned  surface  vehicle  (USV) 
with  a  receiver  was  used  as  a  proxy  for  an  unmanned  underwater  vehicle  (UUV).  The 
results  of  these  field  tests  were  published  in  a  paper  titled,  “Underwater  acoustic 
navigation  using  a  beacon  with  a  spiral  wave  front  [3].” 

To  assess  the  performance  of  the  spiral  wave  front  beacon  in  operational,  ocean 
environments,  it  was  necessary  to  develop  a  technique  to  properly  apply  acoustic 
propagation  models  to  the  novel,  spiral  wave  front  source.  This  was  accomplished  and 
the  techniques  were  published  in  a  paper  titled,  “Acoustic  propagation  from  a  spiral  wave 
front  source  in  an  ocean  environment  [4].” 

FY12-FY13 


The  goal  in  the  last  two  years  of  the  project  was  to  develop  and  test  a  final,  deployable 
version  of  the  spiral  beacon  transducer.  To  move  closer  to  a  fleet-ready  version  of  the 
spiral  beacon,  the  USV  used  in  previous  field  tests  was  to  be  replaced  with  a  UUV  to 
assess  the  performance  of  the  beacon  at  longer  ranges,  various  depths,  and  under  more 
realistic  operating  conditions.  The  modeling  techniques  developed  in  the  previous  year 
were  to  be  implemented  in  acoustic  propagation  codes  and  used  to  assess  the 
performance  of  the  beacon  in  various  operational  environments. 

Accomplishments 

Following  the  field  tests  in  FY  10  and  FY1 1,  Dr.  David  Brown  at  BTech  Acoustics  LLC 
developed  a  final,  production-ready  spiral  beacon  prototype.  To  evaluate  the 
performance  of  this  new  prototype,  a  receiver  was  integrated  into  an  autonomous 


underwater  vehicle  (AUV)  and  this  AUY  was  used  to  evaluate  the  beacon  performance  at 
the  Navy’s  Seneca  Lake  Sonar  Test  Facility  operated  by  NUWC-Newport.  These  tests 
occurred  in  the  summer  of  FY13.  The  spiral  beacon  propagation  modeling  techniques 
were  implemented  using  a  normal  mode  code  and  used  to  evaluate  the  performance  of  the 
beacon  during  the  Seneca  Lake  tests.  The  results  of  the  experiment  and  the  propagation 
modeling  were  presented  at  the  2014  ONR  Unmanned  Maritime  Systems  Technology 
(UMST)  Program  Review  in  Panama  City  Beach,  FL. 

The  spiral  beacon  normal  mode  code  was  also  used  to  evaluate  the  performance  of  the 
spiral  beacon  for  use  in  navigation  of  a  USV-deployed  mine  neutralizer.  The  results  of 
this  assessment,  much  of  which  applies  to  navigation  of  any  UUV  using  the  spiral  wave 
front  beacon,  were  published  in  the  final  report  for  that  project  [5].  While  the  normal 
mode  code  provides  a  high  fidelity  means  of  evaluating  the  spiral  beacon  performance,  it 
is  not  practical  to  use  at  the  center  frequencies  and  bandwidths  for  which  the  final 
prototype  has  been  designed.  Efforts  have  been  underway  to  implement  the  spiral  beacon 
into  the  Navy’s  Sonar  Simulation  Toolset  developed  by  Dr.  Robert  Goddard.  Final 
evaluations  of  this  implementation  have  not  been  completed  at  the  close  of  this  project 
and  will  be  completed  under  a  different  project. 
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A  spiral  wave  front  source  produces  an  acoustic  field  that  has  a  phase  that  is  proportional  to  the  azi¬ 
muthal  angle  about  the  source.  The  concept  of  a  spiral  wave  front  beacon  is  developed  by  combin¬ 
ing  this  source  with  a  reference  source  that  has  a  phase  that  is  constant  with  the  angle.  The  phase 
difference  between  these  sources  contains  information  about  the  receiver’s  azimuthal  angle  relative 
to  the  beacon  and  can  be  used  for  underwater  navigation.  To  produce  the  spiral  wave  front,  two 
sources  are  considered:  a  “physical-spiral”  source,  which  produces  the  appropriate  phase  by  physi¬ 
cally  deforming  the  active  element  of  the  source  into  a  spiral,  and  a  “phased-spiral”  source,  which 
uses  an  array  of  active  elements,  each  driven  with  the  appropriate  phase,  to  produce  the  spiral  wave 
front.  Using  finite  element  techniques,  the  fields  produced  by  these  sources  are  examined  in  the 
context  of  the  spiral  wave  front  beacon,  and  the  advantages  of  each  source  are  discussed. 
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I.  INTRODUCTION 

As  autonomous  underwater  vehicles  become  common 
tools  in  commercial,  military,  and  scientific  applications,  there 
is  an  increasing  need  for  accurate  underwater  navigation  and 
positioning.  This  need  has  led  to  advances  in  the  application 
of  a  wide  range  of  technologies,1’2  yet  acoustic  positioning 
remains  a  central  and  robust  technique  for  underwater  naviga¬ 
tion.  Among  the  acoustic-based  tools  used  for  underwater  nav¬ 
igation  are  long  baseline  systems,  which  acoustically 
determine  the  range  to  two  or  more  transponders  and  then 
determine  the  position  through  triangulation,  and  ultrashort 
baseline  systems,  which  utilize  an  array  of  transducers  to 
determine  both  bearing  and  range  to  a  single  transponder.3 

Recently,  a  new  technique  was  proposed  which  utilizes 
an  acoustic  beacon  that  transmits  two  signals  that  have  a 
phase  difference  that  depends  on  the  bearing  relative  to  the 
beacon.4-6  This  beacon  is  modeled  after  the  very  high  fre¬ 
quency  omnidirectional  range  (VOR)  system  used  in  aircraft 
navigation.7  A  VOR  station  transmits  two  radio-frequency 
signals,  in  the  108-118  MHz  band,  which  an  aircraft  pas¬ 
sively  detects.  The  first  signal,  the  reference  signal,  has  a 
phase  that  is  constant  in  all  directions,  whereas  the  second 
signal,  sent  from  a  second  antenna,  has  a  phase  that  varies 
with  bearing  from  the  station.  The  receiver  on  the  aircraft 
can  then  determine  the  phase  difference  between  the  two  sig¬ 
nals  and,  from  that  difference,  determine  the  bearing  to  the 
station.  The  aircraft  can  stay  on  that  bearing  to  fly  toward  or 
away  from  the  station.  Different  stations  operate  at  different 
frequencies,  which  are  noted  on  aeronautical  charts  available 
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to  the  pilot.  If  there  are  two  VOR  stations  in  range,  the  loca¬ 
tion  of  the  aircraft  can  be  determined  by  triangulation. 

This  technique  has  been  the  primary  means  of  aircraft 
navigation  since  the  1950s  and,  to  the  best  of  our  knowledge, 
an  acoustic  analog  of  this  system  for  underwater  navigation 
has  never  been  developed.  A  likely  reason  for  this  absence  is 
the  difficulty  in  generating  an  acoustic  wave  with  the  required 
angular  phase  dependence.  An  outgoing  acoustic  wave  that 
has  a  phase  independent  of  direction  can  be  commonly  gener¬ 
ated  by  a  uniformly  vibrating  cylindrical  source.  The  problem 
is  in  generating  an  acoustic  wave  with  phase  depending  line¬ 
arly  on  the  angle  such  that  it  changes  by  2n  for  one  circum¬ 
navigation  of  the  source,  and  thus  having  a  continuous,  spiral 
wave  front.  Hence,  we  will  refer  to  this  source  as  a  “spiral 
source”  and  the  cylindrical  source  as  the  “reference  source.” 

The  purpose  of  this  article  is  to  discuss  how  a  spiral  wave 
front  can  be  generated  and  how  a  beacon  composed  of  both  a 
reference  and  a  spiral  source  ideally  would  behave.  Techni¬ 
ques  that  use  the  output  of  these  sources  for  underwater  navi¬ 
gation  will  be  discussed  in  a  separate  paper.  We  begin  in  Sec. 
II  by  examining  two  analytical  examples  of  spiral  and  refer¬ 
ence  sources  in  order  to  develop  a  fundamental  understanding 
of  the  field  generated  by  a  spiral  source.  These  examples  also 
provide  a  starting  point  from  which  the  field  produced  by  a 
spiral  wave  front  beacon  can  be  explored.  In  Sec.  Ill  two 
source  variations  capable  of  producing  spiral  wave  fronts  are 
modeled  using  finite  element,  as  well  as  approximate  analyti¬ 
cal  techniques.  Each  of  these  sources  has  different  advantages 
and  disadvantages,  which  are  discussed  in  Sec.  IV. 

II.  SIMPLE  SPIRAL  SOURCE  AND  BEACON 

To  understand  the  field  produced  by  a  spiral  wave  front 
source  and  how  a  spiral  beacon  can  be  constructed,  we  begin 
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with  a  simple,  analytic  example.  The  simplest  solution  to  the 
three-dimensional  Helmholtz  equation  capable  of  producing 
spiral  wave  fronts  can  be  expressed  in  the  far-held  as 

JkR 

P(R ,  0,  4>)  =  pckQs—cos(p9)  (1) 

K 

where  k  is  the  wave  number  in  the  medium,  p  is  the  density 
of  the  medium,  c  is  the  sound  speed,  Qs  is  the  source 
strength,  and  R ,  0,  and  (j)  are  the  spherical  coordinates 
defined  in  Fig.  1.  This  solution  assumes  that  the  sound  speed 
in  the  surrounding  medium  is  isotropic.  When  p  =  0,  this  so¬ 
lution  reduces  to  a  simple  point  source,  which  will  be  used 
as  a  reference  source.  When  p  =  1,  the  phase  of  the  outgoing 
held  depends  on  cp  and  the  wave  fronts  at  z  =  0  have  a  con¬ 
tinuous,  spiral  structure.  As  0  approaches  ±n/2,  the  phase 
becomes  indeterminate  and,  as  a  result,  there  is  a  null  along 
the  z-axis.  The  structure  of  the  held  propagating  close  to  the 
z-axis  therefore  has  a  structure  similar  to  a  helicoidal  beam.8 

Using  Eq.  (1),  a  spiral  beacon  can  be  modeled  by  com¬ 
bining  the  reference  source  ( p  =  0)  and  the  spiral  source 
(p  =  1).  Ideally,  these  sources  would  be  collocated  in  space 
and  the  phase  difference  between  outgoing  wave  fronts 
would  be  independent  of  0.  If  the  sources  are  hred  simulta¬ 
neously,  the  magnitude  of  the  pressure  held  goes  to  zero  at 
0  =  0  and  </>  =  0  [note  that  a  phase  shift  of  n  was  added  to 
Eq.  (1)  to  place  the  null  at  </>  =  0].  The  normalized  pressure 
magnitude  is  shown  in  the  upper  panel  of  Fig.  2.  This  null  is 
robust  as  a  function  of  0.  At  </>  =  0,  moving  away  from 
0  =  0  to  0  =  ±26°,  the  held  initially  increases  from  zero 
very  slowly,  reaching  only  1%  of  the  held  maximum. 

Although  having  the  spiral  and  reference  source  collo¬ 
cated  may  be  the  ideal  situation,  in  many  practical  beacon 
designs  the  two  sources  cannot  be  collocated  and  likely  will 
be  separated  vertically.  This  situation  is  shown  in  Fig.  1 


FIG.  1.  Spiral  beacon  configuration  and  coordinate  system.  The  spiral 
source  is  denoted  by  Ss  and  the  reference  source  is  denoted  by  Sc.  The  inset 
shows  an  arrangement  of  three  sources:  Sc i  and  Sc 2  are  reference  sources 
and  Ss  is  the  spiral  source. 


where  the  sources  are  separated  by  2Az.  For  the  reference 
source,  again  Eq.  (1)  with  p  =  0  can  be  used  to  describe  the 
pressure  held.  However,  to  account  for  the  displacement 
from  the  xy  plane,  both  0  and  R  are  replaced  by  0C  and  Rc , 
which  are  determined  from  Fig.  1, 

“,«y  = ® 


and 


Rc=R 


cos (0) 
cos(0c) ' 


(3) 


Fikewise,  for  the  spiral  source,  Eq.  (1)  can  be  used  with 
p  =  1  and  0  and  R  replaced  by  9S  and  Rs ,  which  are  deter¬ 
mined  from 


tan  (9S)  =  tan(0)  + 


A  z 

Rcos(9 ) 


(4) 


Az/X=  0 
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FIG.  2.  The  normalized  magnitude  of  the  sum  of  the  pressures  from  the 
simple  source  and  the  spiral  source  each  calculated  using  Eq.  (1)  with  n  =  0 
and  n=  1,  respectively.  The  pressures  were  calculated  at  R  =  10  m  as  the 
separation  between  the  sources  is  increased  from  Az  =  0  to  Az  =  2/3.  The 
white  lines  superimposed  on  the  plots  correspond  to  the  approximate  loca¬ 
tion  of  the  nulls,  as  given  by  Eq.  (7). 
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and 


RS=R 


cos (0) 
cos(^) ' 


(5) 


Unlike  the  beacon  with  the  collocated  sources,  the  dif¬ 
ference  between  the  phases  of  the  two  sources  is  no  longer 
independent  of  0.  For  small  values  of  0,  close  to  the  xy  plane 
and  in  the  far-held,  the  phase  difference  between  the  refer¬ 
ence  source  and  the  spiral  source  can  be  written  linearly  in  0, 

A  z 

AO  =  (f)  —  n  —  An  — -  0.  (6) 

A 

If  the  phase  difference  between  the  reference  and  the  spiral 
source  were  being  used  to  determine  the  bearing  to  the  bea¬ 
con,  this  expression  shows  that  it  is  necessary  to  know  0  or, 
alternatively,  the  depth  and  range  relative  to  the  beacon. 
Without  that  knowledge,  the  expression  given  in  Eq.  (6)  can 
provide  a  means  of  estimating  the  error  in  the  bearing  deter¬ 
mination  due  to  uncertainty  in  the  receiver  depth. 

This  phase  difference  affects  the  beacon  output  when 
the  sources  are  fired  simultaneously,  as  can  be  seen  in  Fig.  2. 
The  difference  in  the  phases  as  a  function  of  0  has  the  effect 
of  “shearing”  the  collocated  pressure  field.  Using  the  phase 
difference  in  Eq.  (6),  the  approximate  location  of  the  null  as 
a  function  of  0  can  be  determined  from 


A  z 

4>n  =  4n-^-6  -  n2n, 


(7) 


This  increases  the  complexity  of  the  beacon  and  the  signal 
processing  at  the  receiver,  but  still  requires  only  a  single 
beacon. 

While  these  models  for  the  sources  may  seem  overly 
simplistic,  they  capture  many  of  the  important  aspects  of  the 
spiral  field  and  the  beacon  performance.  To  see  this,  we  will 
consider  an  analytical  model,  that  is,  closer  in  geometry  to 
the  spiral  sources  and  beacons  that  will  be  discussed  in  the 
next  section. 

Each  of  the  sources  is  embedded  in  an  infinitely  long, 
cylindrical  baffle  and  the  active  area  of  the  source  wraps 
completely  around  the  cylinder  with  a  width  of  h  in  the 
z-direction  as  illustrated  in  Fig.  3.  The  active  area  of  each 
source  vibrates  radially,  but  the  phase  of  the  vibration  is  con¬ 
stant  for  the  reference  source  and  is  proportional  to  </>  for  the 
spiral  source.  To  determine  the  fields  generated  by  these 
sources  we  will  use  the  general  expression  for  the  far-held 
radiation  from  a  rectangular  piston  on  an  infinite  cylinder,9 


p(R,  9,  0) 


PC  e'kR  y'  f  kz) 

n  R  ^  cos  9H'  (ka  cos  9)  ’ 

n=— oo  n\  / 


(12) 


where  Wn  is  the  Fourier  transform  of  the  radial  velocity  of 
the  piston,  kz  —  k  sin  0,  a  is  the  radius  of  the  cylinder,  and  H'n 
is  the  derivative  with  respect  to  the  argument  of  the  first 
Hankel  function  of  order  n. 

If  a  rectangular  piston  has  a  height  h ,  an  angular  width 
2a,  and  a  radial  velocity  amplitude  U  =  Uo  exp(— i/up),  the 
Fourier  transform  yields 


where  n  m  0,  ±1,  ±2, ....  This  result  is  plotted  in  Fig.  2 
where  it  tracks  the  held  minimum  well  for  |0|  <  20°. 

The  conhguration  shown  in  Fig.  1  is  only  one  of  many 
possible  arrangements  of  sources  that  can  be  used  to  con¬ 
struct  a  beacon.  For  this  conhguration,  as  Eq.  (6)  shows, 
once  the  receiver  moves  out  of  the  plane  of  the  beacon,  the 
phase  difference  between  the  sources  becomes  a  function  of 
both  </>  and  0.  This  dependence  can  be  eliminated  by  the 
addition  of  a  second  reference  source  as  shown  in  the  inset 
of  Fig.  1.  In  this  arrangement,  the  phase  difference  between 
the  upper  reference  source,  Sc i ,  and  the  spiral  source  is 

A  z 

A®i  =  </>  —  7i  —  2n— -0,  (8) 

A 

whereas  the  phase  difference  between  the  lower  reference 
source,  Sc 2,  and  the  spiral  source  is 

A  z 

A®2  =  </>  —  71  +  27C— -0.  (9) 

A 

The  bearing  can  now  be  determined  from 

0  =  -  (A®1  +  A®2)  +  71,  (10) 

while  the  grazing  angle  to  the  beacon  can  be  found  from 


U0Gch  sin[(iu  +  n)  a]  sin(£z/z/2)  _ 

W nn\ kz)  —  /  \  77/0'  (1-^) 

n  (ii  +  n)  a  kzh/2 

For  the  sources  depicted  in  Fig.  3,  the  angular  width  of  the 
piston  is  2n  and  a  =  n.  In  Eq.  (13),  WnfJL( a,  kz)  is  zero  for  all 
n  except  when  n  =  —/i,  in  which  case, 


Z 


,  A<D2  -  AO  I 
4tlAz 


(11) 


FIG.  3.  Spiral  beacon  composed  of  two  circumferential  sources  imbedded 
in  an  infinite,  cylindrical  baffle.  The  spiral  source  is  denoted  by  Ss  and  the 
cylindrical  source  is  denoted  by  Sc. 
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W-w( 71,  kz)  =  2£/0 


sm{kzh/2) 

kz 


(14) 


The  sources  are  moved  from  z  =  0  to  z  =  ±Az  by  substitut¬ 
ing  the  relationships  in  Eqs.  (2)-(5).  For  the  spiral  source  (Ss 
in  Fig.  3),  /i  =  1  and  the  only  nonvanishing  term  in  Eq.  (12) 
is  n  =  —  1.  The  far-held  pressure  from  the  spiral  source 
becomes 


terference,  Eq.  (7),  is  superimposed  on  the  cylindrical  source 
interference  in  Fig.  5,  we  see  that  it  tracks  the  interference 
pattern  very  well.  The  equations  for  the  radiation  from  the 
cylindrical  sources  in  Eqs.  (15)  and  (16)  can  be  expanded 
about  0  =  0  and  in  this  small  angle  approximation  it  is 
straightforward  to  show  that  the  phase  difference  between 
the  sources  is  described  by  Eq.  (6).  Likewise,  the  results  for 
the  three  source  beacon  also  apply  to  the  cylindrical  beacon. 


ps(R,  9,  <f>) 


2'U  e‘kRs  s* n(kzsh/2)  e 

°  Rs  nk2S  cos  9S  //«'(&,  cos  0S)  ’ 


(15) 


where  kzs  =  ksin6s.  For  the  reference  source  (Sc  in  Fig.  3), 
H  =  0  and  the  only  nonvanishing  term  is  n  =  0.  The  far-held 
pressure  from  the  reference  source  is  then 


gikRc 

Pc(R,  0,  4>)  =  'ZiUcpc—— 

*<C 


sin(&zc/z/2)  1 

7ikzc cos 9C  //(!)' (ka cos 9C)’ 


(16) 


where  kzc  =  k  sin  0C  and 


fCM 

H[l)\ka) 


Ur  =  -i 


Uo 


(17) 


is  the  radial  velocity  of  the  SC9  dehned  such  that  \pc\  =  \ps\ 
when  0  =  0 

An  example  of  the  phase  of  the  cylindrical  spiral  source  is 
shown  in  Fig.  4.  For  this  example,  R  =»  10  m,  a  m  13  cm, 
/  =  79.5  kHz,  and  h  =  X  =  1.9  cm.  Unlike  the  simple  sources, 
the  cylindrical  sources  each  have  a  main  lobe  which,  for  this 
particular  choice  of  h ,  has  a  vertical  beam  width  (3  dB  down) 
of  58°.  As  a  result,  the  phase  of  the  reference  source  is  no  lon¬ 
ger  constant  over  all  angles,  instead  varying  with  increasing  0. 
The  spiral  source  has  the  same  beam  pattern,  whereas  the  phase 
dependence  in  0  is  superimposed  on  the  phase  ramp  in  </>. 

In  the  beacon  configuration  shown  in  Fig.  3,  when 
Az  =  0  and  the  sources  are  collocated,  the  sum  of  the  fields 
from  the  two  sources  again  produces  a  null  at  the  origin.  As 
Az  increases,  as  in  Fig.  5,  the  interference  pattern  is 
“sheared”  as  it  was  for  the  simple  sources  in  Fig.  2.  In  fact, 
if  the  equation  of  the  field  minimum  for  the  simple  source  in- 
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FIG.  4.  Phase  of  the  cylindrical  spiral  source  at  R  =  10  m  calculated  using 
Eq.  (15).  For  this  source,  h  —  X,  a  =  7.3  cm,/  =  79.5  kHz. 


III.  SPIRAL  WAVE  FRONT  SOURCES 

We  will  focus  on  two  spiral  sources,  each  of  which  uses 
a  different  approach  to  producing  a  wave  front  with  the 
proper  phase  structure.  The  first  source,  which  will  be 
referred  to  as  the  physical-spiral  source,  creates  the  phase  by 
physically  deforming  the  active  element  of  the  source.  The 
second  source,  referred  to  as  the  phased-spiral  source,  uses 
an  array  of  elements,  each  driven  with  a  different  phase,  to 
produce  the  spiral  wave  front. 

This  approach  to  producing  the  desired  phase  variation 
is  similar  to  that  used  to  demonstrate  helicoidal  beam  gener¬ 
ation  in  Ref.  8.  The  first  transducer  in  that  paper  used  a  brass 
ring,  which  was  cut  and  deformed  to  create  a  helical  surface 
on  which  a  flexible  PVDF  sheet  was  attached.  This  physical 
deformation  corresponded  to  the  desired  phase  structure  and 
thus  produced  a  helicoidal  beam.  The  second  transducer 
used  a  four-element  array  and  drove  each  element  90°  out  of 
phase  with  the  previous  element  producing  a  360°  change  in 
phase  around  the  axis  of  the  beam.  This  also  produced  a  heli¬ 
coidal  wave  front  in  the  far-held. 

A.  Physical-spiral  source 

The  geometry  of  the  physical- spiral  source  is  shown  in  Fig. 
6.  The  active  element  is  a  single  piece  that  has  been  wrapped 
around  the  circumference  of  a  backing  of  height  h.  The  backing 
has  been  cut  such  that  the  radius  is  a  function  of  </>, 
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FIG.  5.  The  normalized  magnitude  of  the  sum  of  the  pressures  from  the  cy¬ 
lindrical  reference  and  cylindrical  spiral  source  at  R  =  10  m  when  the  sepa¬ 
ration  between  the  sources  is  Az  =  2/3.  The  white  lines  superimposed  on 
the  plots  correspond  to  the  approximate  location  of  the  null  calculated  for 
the  simple  sources  and  given  by  Eq.  (7). 
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FIG.  6.  Geometry  of  the  physical-spiral  source.  The  phase  ramp  to  produce 
the  spiral  wave  front  is  created  by  wrapping  a  single  active  element  around 
a  backing  whose  radius  has  been  cut  into  a  spiral. 


where  a®  =  am{n  -{-1/2,  and  this  change  in  the  radius  forms  a 
spiral.  Although  this  produces  the  required  phase  change  in 
the  outgoing  wave,  it  also  results  in  a  discontinuity  in  the 
active  element,  with  the  ends  of  the  element  offset  by 
A  a  =  X  from  one  another  at  </>  =  0. 

To  understand  the  effect  of  this  discontinuity  in  the 
active  element,  this  source  was  modeled  using  finite  element 
(FE)  methods.  In  the  FE  model,  the  vibrating  element  had  a 
height  h  equal  to  the  height  of  the  backing.  The  top  face,  bot¬ 
tom  face,  and  the  face  on  the  inside  of  the  discontinuity  were 
treated  as  rigid.  The  pressure  field  generated  by  this  source 
was  calculated  inside  a  cylindrical  volume,  which  sur¬ 
rounded  the  source  and  a  set  of  perfectly  matched  layers 
were  used  to  prevent  reflections  from  the  boundaries  of  the 
computational  domain.10  The  far-field  pressure  was  then 
determined  from  the  pressure  and  normal  velocity  of  the 
field  on  the  boundary  of  the  cylindrical  volume  using  the 
Helmholtz  integral.9 

An  example  of  the  far  field  from  the  physical-spiral 
source  is  shown  in  Fig.  7.  In  this  example,  c  =  1500  m/s, 
h  =  X  —  1.9  cm,  f  =  c/s l  =  78.9  kHz,  and  ao  =  8.3  cm. 
Note  the  horizontal  stratification  in  the  magnitude;  these  dif¬ 
fraction  effects  are  due  to  the  horizontal  edges  along  the  top 
and  bottom  surfaces.  Also,  as  one  might  expect,  the  magni¬ 
tude  of  the  pressure  field  is  not  symmetric  about  the  source 
as  it  was  in  the  examples  in  the  previous  section,  but  rather 
has  pronounced  variations  in  the  vicinity  of  the  discontinuity 
at  </>  =  0°.  In  order  to  understand  how  this  discontinuity 
affects  the  outgoing  field,  we  will  employ  a  rather  crude 
model  of  the  source. 

We  begin  by  dividing  the  vibrating  surface  of  the  source 
into  N  segments  each  with  an  angular  width  of  2a,  where 
a  =  7 i/N.  The  distance  from  the  center  of  the  source  to  the 
center  of  element  /  is  <2/  =  ao  +  (/  —  (N  —  l)/2)  2/A.  We 
will  approximate  each  of  these  elements  as  a  piston  of  height 
h  and  width  2a  on  a  infinite,  cylindrical  baffle  of  radius  at. 
Note  that  although  this  fails  to  capture  the  finite  height  of  the 
FE  model  of  the  source,  the  primary  goal  of  this  approxima¬ 
tion  is  to  capture  the  behavior  of  the  outgoing  field  close  to 
the  xy  plane.  Using  the  expression  for  a  piston  on  a  cylindrical 
baffle,  Eq.  (12),  the  pressure  field  due  to  the  Zth  element  is 
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FIG.  7.  Pressure  (upper)  and  phase  (lower)  calculated  at  R  =  1  m  from  a  FE 
simulation  of  the  physical-spiral  source  shown  in  Fig.  6.  The  pressure  has 
been  normalized  by  the  maximum  pressure  determined  from  the  FE 
calculation. 


(i9) 

R  cos  6H'n [kai cos  0) 


where  </>/  =  </>  —  n  —  (/  +  1/2)  2n/N.  The  total  field  is  then 
approximated  as  the  sum  of  these  individual  elements, 


pt{R,  e,  <f>)  = 


pceikR  in<*-n)Wn o(«,  kz)Fn(k,  9) 

n  R  ^  1  ’  cos0 

n=— oo 


where 


Fn(k,  9) 


^4exp[-m(/  +  I)^] 
H'n(kai  cos  9) 


(20) 


(21) 


As  N  — >  oo,  this  expression  captures  the  ramp  in  the  radius 
as  a  function  of  </>.  At  the  discontinuity,  although  this 
approximation  captures  the  abrupt  change  in  the  radius  at 
that  point,  it  does  not  account  for  the  rigid  face  that  runs 
along  the  x  axis  from  a  =  ao  —  2/2  to  a  =  ao  +  2/2.  In  this 
approximation,  as  the  discontinuity  is  approached  from  ei¬ 
ther  positive  or  negative  </>,  the  border  at  </>  =  0  is  a  cylindri¬ 
cal  baffle,  not  the  rigid  face  or  edge  as  in  Fig.  6. 

The  field  produced  by  this  approximation  to  the  source 
response  is  shown  in  Fig.  8  and,  despite  the  crudeness  of  the 
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FIG.  8.  Pressure  (upper)  and  phase  (lower)  calculated  at  R  =  1  m  by 
approximating  the  field  from  the  physical- spiral  source  by  a  sum  of  pistons 
embedded  in  cylindrical  baffles  using  Eq.  (20).  For  these  results,  N  =  90. 
The  pressure  has  been  normalized  by  the  maximum  pressure  determined 
from  the  FE  calculation  in  Fig.  7. 


model,  it  does  capture  aspects  of  the  overall  structure  of  the 
field  determined  from  the  FE  calculation.  Specifically,  it 
shows  very  similar  variations  in  the  magnitude  near  the 
angular  location  of  the  discontinuity.  As  the  model  does  not 
capture  the  affect  of  the  rigid  face  at  the  discontinuity,  these 
similarities  arise  from  the  difference  in  distance  between  the 
elements  on  either  side  of  the  discontinuity.  For  the  portion 
of  the  element  that  lies  along  (j)  >  0°,  the  radiated  field  must 
travel  an  extra  wavelength  relative  to  the  field  produced  by 
the  element  along  (j)  <  0°,  and  this  extra  distance  results  in 
additional  spreading  loss.  Unlike  the  previous  simulation, 
this  method  exhibits  no  edge  diffraction  effects  due  to  the 
lack  of  top  and  bottom  surfaces. 

To  distinguish  the  effect  of  the  change  in  element  dis¬ 
tance  from  that  of  the  rigid  wall  at  the  discontinuity,  the  field 
along  9  =  0°  is  shown  in  Fig.  9.  For  the  amplitude  compari¬ 
son,  the  fields  compare  very  well  within  <fi  =  ±20°  of  the 
discontinuity  and  both  show  oscillations  in  the  field  beyond 
this  interval  out  to  </>  =  ±90°.  The  oscillations  are  larger  for 
the  approximation  and  are  not  quite  in  phase  with  the  simu¬ 
lation.  This  difference  is  likely  due  to  the  absence  of  the 
rigid  face  at  </>  =  0°  in  the  approximation. 

These  differences  are  more  pronounced  in  the  plot  of 
the  phase  variation,  (5®  =  ®  +  </>  —  n,  where  ®  is  the  phase 
of  the  physical-spiral  output,  shown  in  the  lower  panel  of 
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FIG.  9.  Output  of  the  physical-spiral  source  at  6  =  0°  determined  from  the 
FE  calculation  (solid  line)  and  the  approximation  (dashed  line)  given  by  Eq. 
(20).  The  upper  plot  shows  the  amplitude  of  the  pressure  field  normalized  by 
the  maximum  value  of  the  FE  calculation.  The  lower  plot  is 
<5cp  =  cD  +  (/)  —  7i.  Note  that  the  phase  of  an  ideal  spiral  source  would  be 
<D  =  —(</>  —  7i )  and  hence  <D  +  (j)  —  n  would  be  zero. 


Fig.  9.  For  an  ideal  spiral  source,  ®  =  —(</>  —  n)  and  the 
phase  variation  would  be  zero.  Again  the  oscillations  are 
larger  for  the  approximation  than  the  simulation.  Although 
the  oscillations  in  the  phase  are  due  to  the  change  in  distance 
from  the  elements  on  either  side  of  the  discontinuity,  the 
rigid  face  in  the  simulation  produces  an  overall  increase  in 
phase  for  </>  <  0°  and  overall  decrease  in  phase  for  </>  >  0°. 
In  both  the  simulation  and  the  approximation,  the  phase  dif¬ 
ference  varies  between  ±10°  whereas  the  standard  deviation 
is  3.65°  for  the  approximation,  due  largely  to  the  oscilla¬ 
tions,  and  4.42°  for  the  simulations. 

Because  the  approximation  uses  pistons  on  an  infinitely 
long  cylindrical  baffle,  there  are  pronounced  differences 
between  the  approximation  and  the  simulation  away  from 
the  xy  plane.  In  the  pressure  magnitude  plots,  the  finite  size 
of  the  source  gives  rise  to  oscillations  as  \9\  increases.  These 
are  absent  from  the  approximation  where  the  field  decreases 
smoothly  to  zero  as  \()\  increases.  In  the  phase  plots,  the 
approximation  shows  a  pair  of  phase  dislocations  at 
9  =  ±45°  and  0  =  0°.  Although  difficult  to  see  in  the  gray¬ 
scale  plot  in  Fig.  7,  the  phase  determined  from  the  FE  calcu¬ 
lation  also  exhibits  phase  dislocations,  but  at  higher  grazing 
angles,  9  «  ±60°. 

A  simple  beacon  can  be  constructed  from  the  physical- 
spiral  source  and  a  cylindrical  reference  source.  The  refer¬ 
ence  source  is  chosen  to  be  a  cylindrical  source  similar  in 
construction  to  the  physical-spiral  source  shown  in  Fig.  6, 
but  with  a  constant  radius,  and  thus  no  discontinuity.  The 
spiral  and  reference  are  arranged  with  a  common  axis  in  the 
z-direction  with  the  center  planes  of  the  sources  separated  by 
a  distance  2Az  similar  to  the  diagram  in  Fig.  3. 

As  in  the  simple  source  examples  of  the  previous  sec¬ 
tion,  the  phase  difference  between  the  reference  and  spiral 
sources  can  be  approximated  by  Eq.  (6)  near  0  =  0°.  This 
expression  indicates  that  the  closer  the  sources  are  to  one 
another,  the  dependence  of  the  phase  difference  on  9  is 
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FIG.  10.  The  phase  difference  between  the  beacon  sources  when  Az  =  1.60 
cm  as  a  function  of  pitch  angle  for  three  different  aspect  angles:  4>  =  180° 
(Top  panel),  </>  =  90°  (Middle  panel),  and  $  =  0°  (Bottom  panel).  The  solid 
line  is  the  phase  difference  determined  from  the  FE  calculations  of  the  bea¬ 
con.  The  dashed  line  is  determined  using  Eq.  (6).  The  dotted  line  is  deter¬ 
mined  using  Eq.  (24). 

weaker.  However,  because  of  the  difference  in  the  shapes  of 
the  physical-spiral  and  the  reference  sources,  a  smaller 
source  separation  leads  to  an  increase  in  the  interaction  of 
the  outgoing  field  with  the  structure  of  the  other  source.  In  fi¬ 
nite  element  simulations  of  the  full  beacon,  when  the  separa¬ 
tion  is  A z  =  1.05  cm,  there  are  pronounced  effects  on  the 
outgoing  field  from  both  the  spiral  and  the  reference  source, 
which  are  reduced  significantly  when  the  separation  is 
increased  to  A z  =  1.6  cm. 

The  approximation  for  the  phase  difference  between  the 
reference  and  the  spiral  sources  can  be  improved  by  consid¬ 
ering  the  effect  of  the  finite  radius  of  the  beacon  sources.  In 
Fig.  1,  the  distances  Rc  and  Rs  are  measured  from  the  point 
sources  on  the  z-axis.  For  the  source  geometry  shown  in  Fig. 
6,  the  distances  can  be  approximated  as  originating  from  the 
vertical  center  of  the  active  element  on  the  outer  edge  of  the 
source.  In  this  case,  the  distances  can  be  written  as 

Rc  =  \J (/?  sin  6  -  A z)2+(fi  cos  9  -  a0)2  (22) 
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FIG.  11.  Geometry  of  the  phased-spiral  source.  The  active  element  is 
wrapped  around  a  cylindrical  surface  and  then  etched  into  N  individual  ele¬ 
ments,  each  of  which  is  driven  with  the  appropriate  phase  to  produce  an  out¬ 
going  spiral  wave. 

and 

Rs  =  yj ( R  sin  6  +  A z)2+(fl  cos  6  -  a((/>))2,  (23) 

where  a(</>)  is  given  in  Eq.  (18).  In  the  far- field,  where 
Rsin9^>Az  and  R  cos  9  a((j)),  the  phase  difference, 

AO  =  k{Rc  —  Rs),  can  be  approximated  as 

4-71 

AO  =  ((/)  —  n)  cos  9 - —  Az  sin  9.  (24) 

A 

The  unwrapped  phase  difference  for  the  spiral  beacon  as 
a  function  of  pitch  angle  is  shown  in  Fig.  10  for  three  differ¬ 
ent  aspect  angles  and  compared  to  the  small-angle  approxi¬ 
mation  given  by  Eq.  (6)  and  the  physical-spiral 
approximation  given  by  Eq.  (24).  In  these  comparisons, 
again,  c  =  1500  m/s,  h  =  2  =  1.9  cm,  /  =  cjX  =  78.9  kHz, 
and  ao  =  8.3  cm.  In  all  three  cases,  the  small  angle  approxi¬ 
mation  and  the  physical-spiral  approximation  capture  the  de¬ 
pendence  on  9  within  9  =  ±20°.  For  </>  =  180°  and 
</>  =  90°,  the  physical- spiral  approximation  captures  the 
trend  in  the  phase  difference  over  9  =  ±80°.  At  </>  =  0°,  the 
angle  at  which  the  offset  in  the  physical-spiral  source  is 
located,  neither  approximation  captures  the  phase  difference 
for  steep  angles,  which  is  expected  due  to  the  influence  of 
the  offset. 

The  results  in  Fig.  10  indicate  that  for  small  angles, 
determination  of  the  aspect  angle  using  the  phase  difference 
between  the  spiral  beacon  sources  is  quite  robust  at  any  as¬ 
pect  angle.  For  aspect  angles  opposite  the  spiral  offset  (the 
back  side  of  the  beacon),  the  technique  is  robust  over  nearly 
all  pitch  angles. 

B.  Phased-spiral  source 

For  the  phased-spiral  array  geometry  shown  in  Fig.  11, 
the  active  element  is  wrapped  around  a  cylindrical  backing 
and  is  etched  into  N  individual  elements,  each  of  which  is 
driven  independently.  The  half-angular  width  of  each 
element  is  a  =  n/N.  To  produce  the  spiral  wave  front,  each 
element  /  is  driven  with  a  voltage  that  has  a  constant  ampli¬ 
tude,  but  a  phase  that  depends  on  the  angular  position  of  the 
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element,  producing  a  different  radial  velocity  for  the  surface  These  fields  are  then  summed  to  determine  the  total  field 
of  each  element,  produced  by  the  phased-spiral  source, 


Ui  =  Uq  exp 


(25) 


where  /  =  0,  1, ...,  TV  —  1.  This  discrete  phase  ramp  should 
produce  an  approximately  continuous  phase  ramp  in  the  far- 
field  of  the  source  as  can  be  seen  in  the  FE  calculation  results 
in  Fig.  12  for  a  source  with  =  16.  Again  the  FE  simulation 
treats  each  element  as  a  vibrating  surface,  neglecting  the  ma¬ 
terial  and  electrical  response  of  the  individual  elements  or 
the  backing.  In  this  example,  c  =  1500  m/s,  h  =  2  =  1.9  cm, 
/  =  c// 1  =  78.9  kHz,  and  a  =  8.3  cm.  Although  this  does 
produce  the  desired  phase  ramp,  both  the  amplitude  and 
phase  of  the  field  have  oscillations  in  the  ^-direction  due  to 
the  summation  of  the  fields  from  the  individual  elements. 

To  examine  these  variations  in  the  field,  we  will  approx¬ 
imate  the  output  of  the  phased-spiral  source  by  a  phased 
array  of  pistons  on  a  cylindrical  baffle.  The  field  produced 
by  an  individual  piston  is  given  by  Eq.  (12)  with  </>  replaced 
by  </>/  =  </>  —  7i  —  l(2n/N),  and  the  Fourier  transform  of  the 
radial  velocity  becomes 

w$  (a,  kz)  =  W„o(cc,  kz)  exp f-i/^ .  (26) 
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FIG.  12.  Pressure  (upper)  and  phase  (lower)  calculated  at  R  =  1  m  from  a 
FE  simulation  of  the  phased-spiral  source  shown  in  Fig.  11.  The  pressure 
has  been  normalized  by  the  maximum  pressure  determined  from  the  FE 
model. 


pce^_  (-/)V«(4>-*)  1^,(0 c,  kz) 

n  R  ^  cos  OH'  (ka  cos  9) 

n=— oo  n\  J 


X 


N- 1 


E 


exp 
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2n 
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(27) 


The  summation  over  /  can  be  written  as  a  sum  of  Kronecker 
deltas, 


N- 1 


Eexp 

1=0 


—i(n  +  1)  / 


2n 

W 


—  N  ^  dn(mN- l)j 


(28) 


where 

&n(mN—  1) 


1  if  n  =  mN  —  1 
0  otherwise. 


(29) 


Substitution  of  Eq.  (28)  into  Eq.  (27)  yields 

=  pce^_  ^  (-/)mAf-1e^-^)e-»(0-^)H/(mAf_1)o(a,  kQ 
P  n  R  mAoo  COS  i  (ka  cos  0) 

(30) 


Note  that  when  N  m  l,  there  is  only  a  single  piston  and 
W(mN- i)o  =  0  for  all  m  except  m  =  1.  The  expression  then 
reduces  to  the  solution  for  the  reference  source  given  by  Eq. 
(16),  with  Uc  =  Uo.  When  the  number  of  pistons  becomes 
very  large,  the  derivative  of  the  Hankel  function  becomes 
very  large  for  all  terms  in  the  sum  except  m  =  0  and  the 
expression  reduces  to  the  spiral  source  result  in  Eq.  (15)  as 
expected. 

An  example  of  the  field  calculated  using  this  approxima¬ 
tion  is  shown  in  Fig.  13  for  N  =  16.  Again,  there  are  differ¬ 
ences  between  the  approximation  and  the  FE  calculation  as 
\0\  increases  due  to  the  presence  of  the  cylindrical  baffle  in 
the  approximation.  However,  close  to  6  =  0°  the  approxima¬ 
tion  and  the  FE  calculation  show  a  very  similar  structure, 
specifically  the  oscillations  in  the  magnitude  and  phase  due 
to  the  fields  from  the  individual  elements.  When  9  =  0°,  the 
two  results  compare  very  well  and  the  approximation  will  be 
used  to  investigate  the  source  performance  near  9  =  0° . 

To  predict  the  phased-spiral  source  performance,  we 
will  consider  the  standard  deviation  of  the  normalized  ampli¬ 
tude  and  of  the  phase  variation.  When  used  as  a  beacon, 
these  standard  deviations  affect  the  error  in  the  determina¬ 
tion  of  the  bearing  of  the  receiver  relative  to  the  beacon. 
This  error  will  also  be  affected  by  the  skewness  of  the  ampli¬ 
tude  and  phase  variation,  but  this  will  have  a  smaller  effect 
on  the  bearing  and  we  will  not  consider  it  here. 

As  stated  earlier,  the  oscillations  in  the  field  along  the 
</> -direction  are  due  to  the  interference  of  the  fields  produced 
by  the  different  elements.  The  interference  pattern  is  related 
to  the  beam  patterns  of  the  elements,  which  depend  on  the 
width  of  the  elements  along  the  circumference,  2oca,  as  well 
as  the  wave  number.  This  interference  pattern  therefore 
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FIG.  13.  Pressure  (upper)  and  phase  (lower)  calculated  at  R  =  1  m  by 
approximating  the  field  from  the  phased-spiral  source  by  a  sum  of  phased 
pistons  in  a  cylindrical  baffle  using  Eq.  (30).  The  pressure  has  been  normal¬ 
ized  by  the  maximum  pressure  determined  from  the  FE  model  in  Fig.  12. 

depends  on  both  N  and  the  ka  of  the  source.  This  is  seen  in 
Fig.  14  where  both  standard  deviations  have  been  plotted  as 
a  function  of  ka  with  TV  =  16.  Also  plotted  for  comparison 
are  the  results  of  FE  calculations  for  phased-spiral  sources 
with  different  values  of  ka.  The  FE  and  approximation  com¬ 
pare  well  over  the  entire  range  of  ka  further  supporting  the 
use  of  this  approximation  when  8  =  0° . 

When  ka  is  small,  both  standard  deviations  become  small 
and  the  held  produced  by  the  array  is  well  approximated  by 
the  cylindrical  source  expression  given  by  Eq.  (15).  To  deter¬ 
mine  what  “small  ka ”  means  in  this  context,  we  consider  how 
the  derivative  of  the  Hankel  function  depends  on  ka  cos  8. 
When  \n\  >  ka  cos  0,  where  n  is  the  index  of  the  Hankel  func¬ 
tion,  1  /H'n  becomes  very  small.  Those  terms  in  Eq.  (30), 
which  satisfy  |  mN  —  1|  >  £z/cos  8  can  be  neglected  from  the 
sum.  This  means  that  when  N  —  1  >  ka  cos  8 ,  all  terms  except 
m  =  0  can  be  neglected  and  the  array  is  well  approximated  by 
the  simple  source.  For  the  array  in  Fig.  14,  this  occurs  when 
ka  <  15,  which  roughly  marks  the  beginning  of  the  decrease 
in  the  standard  deviations  as  ka  decreases. 

As  ka  increases,  the  condition  | mN  —  1|  >  ka  cos  8  con¬ 
tinues  to  influence  the  number  of  terms  that  must  be  included 
in  the  sum  in  Eq.  (30)  in  order  to  calculate  the  field.  In  Fig. 
14,  the  standard  deviation  of  the  normalized  amplitude  is 
well  approximated  by  keeping  only  m  =  — 1 , 1  in  the 
sum  for  ka  <31,  whereas  the  phase  variation  is  well 
approximated  when  ka  <27 .  This  difference  between  the 


FIG.  14.  Standard  deviation  of  the  normalized  amplitude  (upper  plot)  and 
the  phase  variation  (lower  plot)  of  the  outgoing  field  at  R  =  1  m  in  the  plane 
of  a  16-element  phased-spiral  source  calculated  using  the  expression  given 
by  Eq.  (30)  as  a  function  of  ka.  For  each  line  only  the  terms  m  =  — M, . . . ,  M 
were  used,  where  M  =  1  (dash-dot  line),  M  =  2  (dashed  line),  and  M  —  10 
(solid  line).  The  vertical  lines  mark  the  positions  of  the  minima  in  the  stand¬ 
ard  deviations  that  are  plotted  in  Fig.  15.  The  results  determined  from  FE 
calculations  are  shown  as  open  circles. 

two  deviations  is  due  to  the  presence  of  a  decrease  in  phase 
variation,  which  occurs  at  ka  =  28.  This  decrease  means  that 
the  \m\  =2  terms  can  have  an  affect  on  the  phase  variation 
even  though  its  value  is  small.  These  minima  occur  at  differ¬ 
ent  values  of  ka  for  both  standard  deviations  and  in  some 
cases  are  almost  nulls,  meaning  that  the  phase  or  amplitude 
are  nearly  equal  to  that  of  the  simple  spiral  source.  The  field 
itself  never  approaches  that  of  the  simple  spiral  source,  how¬ 
ever,  as  a  local  minimum  for  the  phase  variation  corresponds 
to  a  local  maximum  for  the  amplitude  and  vice  versa.  This 
does  mean,  however,  that  a  signal  processing  technique  that 
only  depends  on  either  the  amplitude  or  the  phase  of  the 


_ i _ i _ i _ i _ 

10  12  14  16  18  20 
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FIG.  15.  Values  of  ka  for  the  local  minima  in  the  standard  deviations  of  the 
normalized  amplitude  (solid  lines  connecting  open  circles)  and  the  phase 
variation  (dashed  lines  connecting  open  diamonds)  as  a  function  of  the  num¬ 
ber  of  elements  in  the  phased-spiral  source.  The  lower  dotted  lines  denotes 
the  condition  N  —  1  =  ka  and  the  upper  dotted  line  denotes  the  condition 
2N  —  1  =  ka. 


3638  J.  Acoust.  Soc.  Am.,  Vol.  129,  No.  6,  June  201 1  B.  T.  Hefner  and  B.  R.  Dzikowicz:  Spiral  wave  front  beacon 


Downloaded  14  Jun  2011  to  128.95.76.64.  Redistribution  subject  to  ASA  license  or  copyright;  see  http://asadl.org/journals/doc/ASALIB-home/info/terms.jsp 


source  field  can  perform  better  if  the  source  is  designed  to 
operate  at  one  of  these  minima. 

The  position  of  these  minima  depends  on  the  number  of 
elements  in  the  array.  If  the  number  of  elements  is 
decreased,  the  minima  shift  toward  lower  values  of  ka.  This 
can  be  seen  in  Fig.  15  where  the  locations  of  the  minima  for 
both  standard  deviations  are  plotted  as  a  function  of  N.  For 
all  values  of  N,  the  minima  are  the  lowest  when  they  occur 
in  the  range  N  <  ka  —  1  <  IN.  The  reason  for  this  can  be 
seen  in  Fig.  14  where  this  range  is  16  <  ka  —  1  <32.  For  ka 
below  this  range,  the  standard  deviations  are  dominated  by 
the  decrease  that  was  discussed  earlier.  Above  this  range,  the 
addition  of  higher  order  terms  in  the  sum  tends  to  fill  in  the 
minima. 

IV.  DISCUSSION 

A  navigational  beacon  can  be  constructed  using  a  source 
capable  of  producing  a  signal  with  phase  that  varies  with  as¬ 
pect  and  a  reference  source  with  phase,  that  is,  constant  with 
aspect.  In  order  to  design  and  optimize  the  performance  of 
such  a  beacon,  the  field  produced  must  be  carefully  under¬ 
stood.  The  modeling  techniques  presented  here  provide  the 
means  of  exploring  different  geometries  and  arrangements 
of  sources. 

The  two  spiral  sources  examined  in  Sec.  Ill  each  use 
very  different  techniques  to  produce  an  approximately  spiral 
wave  front  in  the  far-field.  The  physical- spiral  source  is  the 
simplest  of  the  two  sources  and  requires  only  a  single  ampli¬ 
fier  and  function  generator.  The  operating  frequency  of  this 
source  is  determined  by  the  size  of  the  discontinuity  and  the 
sound  speed  of  the  surrounding  medium.  The  output  of  the 
source  is  strongly  affected  by  the  presence  of  the  discontinu¬ 
ity  at  (j)  =  0°.  This  effect  can  be  minimized  by  either  increas¬ 
ing  the  minimum  radius  of  the  source  or  decreasing  the 
wavelength  of  the  outgoing  field,  and  hence  the  size  of  the 
discontinuity,  at  the  time  of  construction.  Thus  the  ratio, 
A/ao,  must  be  made  as  small  as  possible,  within  the  opera¬ 
tional  requirements  of  the  beacon,  in  order  to  create  the  opti¬ 
mal  spiral  wave  front. 

The  phased-spiral  source  does  not  face  the  same  con¬ 
straints  in  its  operation.  As  each  element  is  driven  individu¬ 
ally,  the  operating  frequency  can  be  changed  while  still 
producing  a  spiral  wave  front  as  long  as  the  phase  difference 
between  the  elements  is  maintained.  This  flexibility,  how¬ 
ever,  comes  at  the  cost  of  increased  complexity  in  the  elec¬ 
tronics  required  to  operate  the  array.  Although  the  number  of 
elements  used  in  the  array  may  be  constrained  by  the  com¬ 
plexity  of  the  electronics,  it  may  be  possible  to  chose  the  fre¬ 
quency  or  radius  such  that  the  operating  ka  corresponds  to 
the  smoothest  phase  ramp,  represented  by  the  minima  shown 
in  Fig.  15.  The  choice  of  minima  will  depend  on  the  signal 
processing  requirements  of  the  beacon  receiver,  either  mini¬ 
mizing  the  phase  variation  or  the  amplitude  variations  of  the 
outgoing  wave  front. 
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Both  sources  can  be  used  with  one  or  two  reference 
sources,  as  discussed  in  Sec.  II,  but  the  phased-spiral  source 
requires  only  a  single  source  to  create  a  beacon.  By  driving 
each  element  with  the  proper  phase,  the  spiral  wave  front 
can  be  generated,  and,  by  driving  each  element  with  the 
same  phase,  the  reference  wave  front  can  also  be  produced. 
This  again  increases  the  complexity  of  the  electronics  but 
produces  a  beacon  for  which  Az  =  0  and  the  phase  differ¬ 
ence  between  the  sources  is  AO  =  (ft  for  all  8. 

Although  the  phased-spiral  source  can  be  optimized  to 
reduce  the  deviation  from  a  true,  spiral  wave  front,  the 
amount  of  deviation  that  is  tolerable  for  underwater  naviga¬ 
tion  remains  to  be  established.  The  errors  calculated  both 
with  the  numerical  and  analytical  models,  as  in  Fig.  14,  rep¬ 
resent  the  minimum  theoretical  error  of  a  real  beacon.  Over¬ 
all  these  errors  are  small,  between  0.1°  and  6°,  for  the 
geometries  and  wavelengths  studied  here.  A  number  of  tech¬ 
niques  to  utilize  the  spiral  beacon  are  currently  being  tested 
to  assess  the  effect  of  the  phase  variation  from  pure  spiral, 
due  to  both  the  environment  and  the  spiral  wave  front  source 
itself,  on  the  navigation  accuracy.11  Future  work  will  also 
have  to  consider  the  effectiveness  of  the  proposed  sources  in 
terms  of  their  transmit  source  level,  operational  frequency, 
bandwidth,  efficiency,  power  factor,  and  beam  patterns. 
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A  spiral  wave  front  source  generates  a  pressure  field  that  has  a  phase  that  depends  linearly  on  the 
azimuthal  angle  at  which  it  is  measured.  This  differs  from  a  point  source  that  has  a  phase  that  is 
constant  with  direction.  The  spiral  wave  front  source  has  been  developed  for  use  in  navigation; 
however,  very  little  work  has  been  done  to  model  this  source  in  an  ocean  environment.  To  this  end, 
the  spiral  wave  front  analogue  of  the  acoustic  point  source  is  developed  and  is  shown  to  be  related 
to  the  point  source  through  a  simple  transformation.  This  makes  it  possible  to  transform  the  point 
source  solution  in  a  particular  ocean  environment  into  the  solution  for  a  spiral  source  in  the  same 
environment.  Applications  of  this  transformation  are  presented  for  a  spiral  source  near  the  ocean 
surface  and  seafloor  as  well  as  for  the  more  general  case  of  propagation  in  a  horizontally  stratified 
waveguide.  © 2012  Acoustical  Society  of  America.  [DOI:  10.1121/1.3682045] 
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I.  INTRODUCTION 

The  navigation  of  a  vehicle  using  a  spiral  wave  front 
beacon  utilizes  the  phase  differences  between  the  spiral 
wave  front  and  reference  sources  that  comprise  the 
beacon.1-3  The  spiral  source  transmits  a  pulse  whose  phase 
is  a  function  of  the  azimuthal  angle  relative  to  the  source 
while  the  reference  source  transmits  a  pulse  whose  phase  is 
constant  with  angle.  By  comparing  the  phases  of  the  two 
pulses,  the  vehicle  can  determine  the  direction  to  the  beacon. 
This  technique  complements  existing  navigation  methods 
such  as  long  baseline  or  ultra-short  baseline,4  with  the  added 
advantage  that  it  can  be  implemented  using  a  single  beacon 
transmitting  to  a  vehicle  with  a  single  receiver. 

Two  types  of  sources  capable  of  producing  the  spiral 
wave  front  have  recently  been  modeled  and  constructed  to 
test  this  navigation  technique.  These  sources  use  two  differ¬ 
ent  methods  to  produce  the  spiral  wave  front.2  The  first,  the 
“physical- spiral”  transducer,  uses  a  single,  active  element 
formed  around  a  spiral  backing.  The  second  design,  the 
“phased- spiral”  array,  uses  an  array  of  elements  each  driven 
with  the  appropriate  phase  to  produce  the  spiral  wave  front. 
Each  of  these  sources  was  modeled  in  the  free  field  to  exam¬ 
ine  how  well  it  produced  a  spiral  wave  front.2 

Recent  experiments  have  emphasized  that  the  perform¬ 
ance  of  the  spiral  beacon  is  affected  not  only  by  the  design 
of  the  source,  but  also  by  the  environment  in  which  it  oper¬ 
ates.5  As  range  from  the  beacon  increases,  it  becomes  neces¬ 
sary  to  understand  how  the  ocean  waveguide  affects  the 
phase  information  carried  by  the  pulses.  Although  there  are  a 
plethora  of  techniques  to  model  the  propagation  of  sound 


from  a  point  source  in  the  ocean,  due  to  the  relatively  recent 
development  of  the  spiral  beacon,  there  has  been  no  work 
modeling  the  propagation  of  sound  from  a  spiral  source.  The 
purpose  of  this  manuscript  is  to  understand  how  the  spiral 
source  and,  by  extension  the  spiral  beacon,  can  be  modeled 
in  the  ocean  environment. 

Understanding  the  propagation  of  sound  from  the  spiral 
source  begins  with  the  development  of  the  spiral  analog  to 
the  ideal  point  source.  The  pressure  field  from  an  ideal  spiral 
source  is  shown  in  Sec.  II  to  be  related  to  an  arrangement  of 
point  sources  and  this  relationship  is  exploited  in  Sec.  Ill  to 
solve  problems  of  reflection  from  the  sea  surface  using  the 
method  of  images  and  from  the  seafloor  using  plane  wave 
decomposition.  In  Sec.  IV,  this  technique  of  transforming 
the  field  of  a  point  source  into  a  spiral  source  is  shown  to  be 
applicable  to  problems  of  propagation  in  a  waveguide  with 
arbitrary  sound  speed  profiles.  Finally  in  Sec.  V,  further 
applications  and  limitations  of  this  approach  are  considered. 

II.  THE  SIMPLE  SPIRAL  SOURCE 

The  acoustic  spiral  wave  front  source  can  be  constructed 
from  two  acoustic  dipoles  driven  90°  out  of  phase  from  one 
another.  The  dipoles  are  located  at  the  origin  and  oriented 
such  that  the  point  sources  that  make  up  the  dipoles  lie  along 
the  x  and  y  axes  as  shown  in  Fig.  1 .  The  pressure  generated 
by  the  first  dipole,  which  lies  along  the  x  axis,  can  be 
expressed  as  the  sum  of  the  field  from  the  two  point  sources,6 


Px(r,e,(j))  =  -ip0ckQs 


G 


,  (1) 
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where  p0  is  the  density,  c  is  the  sound  speed,  k  is  the  wave- 
number,  Qs  is  the  source  strength,  dx  =  dx  is  the  vector  con¬ 
necting  the  two  sources  along  the  x  axis,  and 
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FIG.  1.  Coordinate  system  and  arrangement  of  sources  (open  circles)  used 
to  produce  a  spiral  wave  front. 


G(rV) 


eik\r-r'\ 

|r  -  r'| 


(2) 


which  corresponds  to  the  far-held  source  given  in  Eq.  (1)  of 
Ref.  2  with  Ds  =  Qs/k. 

In  order  to  determine  if  the  expression  given  by  Eq.  (7) 
is  the  simplest  spherical  wave  front  solution  to  the  Helm¬ 
holtz  equation,  the  outgoing  held  can  be  written  in  terms  of  a 
multipole  expansion, 

cx)  n 

p(r,  8,<t>)  =  C™hn\kr)Yn(8>  0)>  (10) 

n= 0  m=—n 

where  hft  are  spherical  Hankel  functions  of  the  hrst  kind 
and  Y ™  are  the  spherical  harmonics.6  Making  use  of  the 
relations 

7;  =  ik2h[l\kr)  (11) 


is  the  Green’s  function  for  a  point  source  located  at  r'.  As 
the  distance  between  the  sources  becomes  vanishingly  small, 
the  pressure  becomes 


Px{r,  9,  <f>)  =  -ip0ckDs 


dG 

dx' 


(3) 


where  Ds  =  Qsd  is  the  dipole  strength.  Expanding  the  partial 
derivative  of  the  Green’s  function, 


8G 

dx' 


r'=0 


8G8R 
8R  dx' 


r'=0 


dG  a  * 
=  — —  cos  0  cos  </>, 

dr 


where  R  =  Ir  —  r'l,  the  pressure  can  be  expressed  as 


px{r,8,$)  =  -ip0ckDs 


-7—  cos  0  cos  </>. 
dr 


(4) 


(5) 


The  dipole  along  the  y  axis  can  be  expressed  in  a  form  simi¬ 
lar  to  Eq.  (3),  with  the  derivative  of  the  Green’s  function 
taken  with  respect  to  /  and  the  expression  multiplied  by  i  to 
account  for  the  phase  shift, 


Py{r,  9,  <f>)  =  p0ckDs 


-7— cos  8  sin  </>. 
dr 


(6) 


The  total  pressure  from  the  sum  of  these  two  dipoles  is 


dG 

ps{r,  9, 4>)  =  -ip^ckD^  —  e1^  cos  9, 


(7) 


which  has  the  angular  phase  dependence  required  for  a  spiral 
wave  front.  In  the  far  field,  the  derivative  of  the  Green’s 
function  is  approximately 


8G 

dr 


and  the  spiral  source  becomes 

ps(r  >  2,  0,  </>)  =  p0ck2Ds- — e1^  cos  0, 


(8) 


(9) 


and 


el<^  cos  0  =  — 


(M), 


(12) 


the  expression  in  Eq.  (7)  can  be  shown  to  correspond  to  the 
n  =  1 ,  m  =  1  multipole  term  with 


(13) 


This  is  the  smallest  term,  and  hence  the  simplest  source,  that 
produces  a  spiral  wave  front.  Note  that  this  is  not  the  only 
multipole  term  that  has  the  spiral  wave  front  phase  depend¬ 
ence.  Any  term  for  which  \m\  =  1  will  have  the  proper  phase 
dependence  since  Y ™  oc  eimc however,  the  0  dependence 
will  be  much  more  complicated  when  n  >  1 . 

For  the  spiral  wave  front  beacons  discussed  in  Ref.  2, 
the  source  strengths  were  chosen  such  that  the  magnitudes  of 
the  fields  were  equal  in  the  far  field  and  when  0  =  0.  If  the 
point  source  is  expressed  as 

eikr 

Po  =  -ip0ckQs  — ,  (14) 

r 


the  far-held  magnitude  of  the  spiral  source  will  be  the  same 
as  that  of  the  point  source  when  the  dipole  strength  is 


Ds  =  -jQs ■  (15) 

k 

With  this  substitution,  the  spiral  source  in  Eq.  (7)  simplifies 
to 


dG 

Ps  =  -p0cQs  -7^7  e1^  cos  9,  (16) 

which  has  the  added  benefit  of  showing  the  explicit  depend¬ 
ence  of  the  held  on  the  source  strength.  The  spiral  source 
can  also  be  written  in  terms  of  the  point  source  by  solving 
Eq.  (14)  for  the  Green’s  function  and  substituting  the  expres¬ 
sion  into  Eq.  (16), 
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(26) 


Ps  =  ~I~^re,<l>  COSd~  (17) 

For  the  point  source,  the  directivity  factor7  is  unity  while  for 
the  spiral  source,  the  directivity  factor  is 

6(0,0)=  cos  0.  (18) 


Ps  =  -i^Ane,(l> 

n 


Q  gikRn 

dp  Rn 


Hence,  any  problem  that  can  be  solved  for  a  point  source 
using  the  method  of  images  or  a  plane  wave  decomposition, 
can  also  be  solved  for  a  spiral  source. 


A  more  complex  source,  such  as  a  cylindrical  transducer  or  a 
small  array,  can  be  approximated  as  the  product  of  a  point 
source  and  the  directivity  factor,  61?  of  the  complex  source, 

P\  =pobi(6,<t>).  (19) 

If  the  directivity  factor  does  not  depend  on  0,  b\  (0,  0)  =  bx 
(0),  then  using  Eqs.  (17)  and  (19),  the  spiral  source  can  be 
written  as 


Ps  =  I  (9)  cos  9,  (20) 

which  has  a  directivity  factor 

b(9)  =  bi(0)  cos  9.  (21) 

This  is  useful  for  modeling  sound  propagation  from  spiral 
sources  such  as  the  physical- spiral  transducer  or  the  phased- 
spiral  array  discussed  in  Ref.  2. 

The  spiral  wave  front  source  given  by  Eq.  (16)  is 
expressed  in  spherical  coordinates.  For  most  problems  in 
ocean  acoustics,  it  will  be  useful  to  rewrite  this  solution  in 
cylindrical  coordinates.  This  is  straightforward  if  we  recog¬ 
nize  that 


dG  dG  dr  dG 

dr  C°S  dr  dp  dp  ’ 


(22) 


where  p  =  r  cos  9  is  the  horizontal  range,  and  express  the  spi¬ 
ral  source  in  cylindrical  coordinates  as 

dG 

Ps{p,z,4>)  =  -p^cQs-^e1^  (23) 

or,  in  terms  of  a  point  source  with  the  same  source  strength, 

Ps(P,  Z,  </>)  =  -  (24) 

This  last  expression  is  the  most  useful  for  solving  problems 
in  ocean  acoustics.  Consider  a  solution  to  the  Helmholtz 
equation  for  a  given  set  of  boundary  conditions  that  can  be 
expressed  as  the  sum  of  point  sources  at  p'  =  0, 


pikRn 

^  =  (25) 

n  n 

where  R2  =  p2  +  (z  —  z^)2  and  An  is  the  pressure  amplitude. 
If  those  point  sources  are  replaced  by  equivalent  spiral  sour¬ 
ces  using  Eq.  (24),  the  solution  to  the  Helmholtz  equation 
for  the  given  boundary  conditions  for  these  spiral  sources 
becomes 


III.  A  SPIRAL  SOURCE  NEAR  AN  INTERFACE 

As  an  example  of  the  application  of  Eq.  (26),  we  will 
begin  by  considering  the  well-known  problem  of  a  point 
source  near  the  ocean  surface.  For  a  point  source  at  a  depth, 
zi,  below  the  surface,  the  held  can  be  expressed  as  the  sum 
of  the  fields  from  the  point  source  and  an  image  source 
located  at  —z\  above  the  surface,8 

gikR  gikR\ 

Po  =  A0  — —  A0  ,  (27) 

where  R  =  p2  +  (z  —  zi)2  and  R\  =  \J  p2  +  (z  +  zi)2.  The 
transmission  loss  in  the  x-z  plane  for  a  point  source  at 
Zi  =  5m  with /=  1.5 kHz  is  shown  in  the  upper  panel  of 
Fig.  2  for  a  medium  with  an  isotropic  sound  speed  of 
cm  1500 m/s.  The  transmission  loss  exhibits  the  familiar 
Lloyd’s  mirror  effect  produced  by  the  interference  of  the 
point  source  and  its  image.9  At  distances  much  greater  than 
the  source  depth,  the  pressure  can  be  approximated  by 

gikRo 

Po(Ro  >  Zi)  «  -2/A0— — sin(foi  sin0o),  (28) 

where  Ro  =  \/ p1  +  z1  is  the  distance  from  the  origin  and 
sin  0O  =  z/^o  is  the  grazing  angle. 

The  field  produced  by  a  spiral  source  near  the 
ocean  surface  can  be  determined  by  applying  Eq.  (24)  to 
Eq.  (27), 


Range  (m) 


FIG.  2.  Transmission  loss  for  a  point  source  (upper  panel)  and  a  spiral 
source  (lower  panel)  below  a  pressure  release  surface.  The  pressure  release 
surface  is  at  the  top  of  each  plot  at  0  m  and  the  sources  have  been  placed  5  m 
below  the  surface.  For  each  source, /=  1.5  kHz  and  c  =  1500  m/s. 
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cos  0 


iAo  , 
-  —  I 

dRi  1 


(29) 


where  cos  0  =  p/R  and  cos  9X  =  p!R\.  The  transmission  loss 
in  the  x-z  plane  for  a  spiral  source  at  z\  =  5  m  is  shown  in 
the  lower  panel  of  Fig.  2  using  the  same  parameters  as  the 
point  source  in  the  upper  panel.  As  expected,  the  most  pro¬ 
nounced  differences  between  the  spiral  and  point  sources  are 
along  the  z  axis  near  p  =  0.  The  cos  9  dependence  of  the  field 
produces  a  null  along  the  z  axis,  leading  to  a  reduction  in  the 
intensity  of  the  interference  fringes  above  and  below  the 
source. 

To  determine  the  pressure  in  the  far  field,  the  approxi¬ 
mations  used  to  determine  the  far-field  approximation  for  the 
point  source  can  be  applied  to  Eq.  (29)  or  Eq.  (24)  can  be 
applied  directly  to  Eq.  (28).  In  both  cases,  the  far-field  pres¬ 
sure  for  the  spiral  source  becomes 


ps(Ro  »  zi) 


—  2iAoe‘^ 


eikR0 

- [sin  (kz\  sin  0O) 

Ro 


i  —  cos(£zi  sin  0q)  sin  0O 
Ro 


cos  0o,  (30) 


10 

Range  (m) 


1000 


FIG.  3.  (a)  Transmission  loss  at  the  depth  of  a  point  source  (solid  line)  and 
a  spiral  source  (dashed  line)  located  5  m  below  a  pressure  release  surface, 
(b)  The  phase  difference  between  a  point  source  and  a  spiral  source  meas¬ 
ured  at  the  source  depth.  The  dashed-dotted  line  is  the  phase  difference 
given  by  Eq.  (36). 


For  small  grazing  angles  and  with  the  source  located  very 
near  the  surface,  sin  0O  «  0O, tan  (kzx  sin  0O)  ~  kzi90,  and  the 
phase  becomes 


For  both  the  spiral  and  the  point  sources,  the  far-field  pres¬ 
sure  is  maximal  when  the  condition 


=  kRo  +  0  —  tan  1  (kRo). 


(34) 


kz\  sin 0O  =  (A  +  1/2)ti,  A  =  0, 1,2, ...  (31) 

is  satisfied.  For  the  far  field  of  the  point  source,  the  acoustic 
pressure  is  zero  when 


Taking  the  phase  of  the  far-field  pressure  for  the  point  source 
from  Eq.  (28), 

$0  =  kR0  -  tc/2,  (35) 


kz\  sin  0o  =  A 7i. 


(32)  the  phase  difference  becomes 


The  nulls  in  the  field  of  the  spiral  source  are  shifted  slightly 
in  angle  relative  to  the  point  source  nulls  and  described  to  a 
good  approximation  by  the  condition  given  in  Eq.  (32).  This 
can  be  seen  clearly  for  ranges  greater  than  20  m  in  Fig.  3(a), 
where  the  transmission  loss  at  the  source  depth  is  plotted  as 
a  function  of  range. 

Note  that  in  Fig.  3(a),  the  near-field  pressure  of  the  spi¬ 
ral  source  diverges  significantly  from  that  of  the  point 
source.  For  the  point  source,  the  near-field  pressure  goes  as 
l/R.  For  the  spiral  source,  the  partial  derivative  of  the 
Green’s  function  leads  to  a  l/R2  dependence  in  the  near  field. 
The  pressure  in  the  far  field  for  the  spiral  source  was  chosen 
to  match  that  of  the  point  source  and,  as  a  consequence,  the 
near-field  pressure  of  the  spiral  source  is  larger  than  that  of 
the  point  source. 

With  these  results,  it  is  possible  to  model  the  perform¬ 
ance  of  a  spiral  wave  front  beacon  operating  near  the  ocean 
surface.  A  spiral  wave  front  beacon  consists  of  a  reference 
source  (point  source)  mounted  2Az  above  the  spiral  wave 
front  source  (see  Fig.  1  in  Ref.  2).  The  phase  of  far-field 
pressure  from  the  spiral  source  is 


^  tan(£Ro  +  0)zi  sin 0o  —  Ro  tan(£zi  sin0o) 

5  z\  sin0o+Ro  tan(£Ro  +  0)tan(£zisin0o)  ’ 


(33) 


AO  =  0  —  tan  l(kRo)  +  7t/2 

Ki4>  +  Tp{/~\e  O’  (36) 

which  depends  very  weakly  on  both  the  grazing  angle  and 
range  when  kp  1.  This  differs  from  the  phase  difference 
for  the  spiral  beacon  operating  in  free  space,2 

A  z 

AO  =  0-  47 1  —  0,  (37) 

A 

where  X  is  the  wavelength  of  the  acoustic  field  generated  by 
the  sources  and  0  «  0O  under  the  approximations  used  in 
Eq.  (34). 

The  far-field  phase  difference  between  the  collocated 
reference  and  spiral  sources  given  by  Eq.  (36)  is  compared 
in  Fig.  3(b)  to  the  exact  phase  difference  determined  from 
Eqs.  (27)  and  (29).  The  far-field  phase  difference  follows  the 
overall  trend  in  the  exact  solution  over  all  ranges.  At  ranges 
less  than  100  m,  the  exact  solution  has  pronounced  oscilla¬ 
tions  that  correspond  to  the  minima  in  the  transmission  loss. 
At  these  points,  the  direct  and  reflected  arrivals  destructively 
interfere  and  these  minima  become  nulls  as  the  range 
increases.  The  discrepancies  in  the  phase  of  the  reference 
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and  spiral  sources  become  more  pronounced  at  these  points 
and  this  produces  the  oscillations  in  the  phase  difference. 

It  may  seem  that  this  kind  of  interference  effect  might 
make  it  difficult  to  use  the  phase  difference  for  navigation. 
In  practical  applications,  the  reference  and  spiral  sources 
send  out  short,  broadband  pulses  and  not  the  continuous 
waves  that  are  discussed  here.  At  very  short  ranges,  the 
direct  and  reflected  arrivals  can  be  distinguished  from  one 
another.  The  interference  becomes  an  issue  at  longer  ranges 
where  the  difference  in  travel  times  between  the  direct  and 
reflected  paths  is  less  than  the  pulse  length.  For  these  ranges, 
the  mean  phase  difference  over  the  frequency  band  of  the 
pulse  can  be  used  with  the  result  being  approximately  the 
phase  difference  given  by  Eq.  (36). 

In  order  to  study  the  performance  of  a  spiral  source  near 
the  seafloor,  we  need  to  know  the  plane  wave  expansion  of 
the  spiral  wave  front  source.  The  expansion  of  a  spherical 
wave  into  plane  waves  is  well  known  and  can  be  written  in 
Cartesian  coordinates  as8 
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where  kz  =  (k2  —  k2  —  k2)  '  and  the  plus  (minus)  sign 
before  kzz  corresponds  to  plane  waves  traveling  in  the  posi¬ 
tive  (negative)  z-direction.  In  cylindrical  coordinates,  the 
plane  wave  expansion  becomes 


(42) 


where  V(kz )  is  the  reflection  coefficient,  kz  =  kcos9,  and  6 
grazing  angle.  When  the  reflection  coefficient  is  independent 
of  angle,  V  can  be  brought  outside  of  the  integral  and  the  so¬ 
lution  reduces  to  that  of  an  image  spiral  source  located  at 
z  =  — z0  multiplied  by  the  reflection  coefficient. 

The  more  complicated  case  of  a  reflection  coefficient 
with  angular  dependence  can  be  solved  using  the  same  tech¬ 
niques  and  approximations  that  are  used  to  solve  for  the 
reflection  of  a  point  source.9  For  example,  to  solve  the  field  at 
distances  large  compared  to  the  wavelength,  the  integral  can 
be  solved  using  the  stationary  phase  approximation  to  yield 
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where  R\  is  the  distance  from  the  receiving  point  to  the 
image  source  location. 

IV.  A  SPIRAL  SOURCE  IN  A  WAVEGUIDE 

The  plane  wave  expansion  of  the  spiral  wave  front 
source  can  also  be  applied  to  propagation  in  a  waveguide 
and,  by  extension,  the  expression  given  by  Eq.  (24)  can  be 
used  to  transform  waveguide  solutions  for  a  point  source 
into  solutions  for  a  spiral  source. 
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As  an  example,  consider  a  point  source  in  a  shallow 
water  waveguide  with  a  flat,  pressure  release  surface  and  a 

(39b) 

flat,  rigid  bottom.  The  solution  for  a  point  source  can  be 
expressed  as  the  sum  of  the  field  from  the  source  and  an  infi¬ 
nite  number  of  image  sources  as 

where  £  =  x  k1  +  k%,  /0  (x)  is  the  Bessel  function  of  the  first 


kind,  and  H ^  (x)  is  the  Hankel  function  of  the  first  kind. 

Using  these  expressions,  it  is  straightforward  to  deter¬ 
mine  the  plane  wave  expansion  of  the  spiral  source.  Apply¬ 
ing  Eq.  (24)  to  Eq.  (39)  yields 
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where  Bn=Jn  or  IlJ  .  was  used.  In  Eq.  (40a),  the  Bessel 
function  of  first  order  ensures  that,  since  the  phase  is  indeter¬ 
minate  at  p  =  0,  the  field  will  have  a  null  along  the  z  axis.10 

Using  the  plane  wave  expansion  of  the  spiral  source,  it 
is  possible  to  determine  the  field  scattered  from  a  flat  fluid  or 
elastic  interface  by  the  spiral  source.  For  a  source  located  at 
z0  above  the  interface,  the  reflected  field  becomes 


z0  is  the  source  depth,  and  h  is  the  depth  of  the  waveguide. 
Using  Eq.  (24),  the  spiral  source  solution  becomes 
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where  cos  6nj  =  p/Rnj  and  j  =  1,  2,  3,  4. 

The  solution  given  by  Eq.  (44)  can  be  expressed  in 
terms  of  normal  modes  by  first  expressing  the  source  and 
image  sources  in  terms  of  their  plane  wave  expansions.  After 
some  manipulations,  the  details  of  which  can  be  found  in 
Chapter  5  of  Ref.  8,  the  field  can  be  expressed  as  a  pair  of 
integrals.  Evaluation  of  those  integrals  reduces  the  point 
source  solution  to  a  sum  of  normal  modes, 
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and  =  \fk2  —  o%.  A  similar  analysis  can  be  used  for  the 
spiral  source  solution  in  Eq.  (49),  but  the  normal  mode  solu¬ 
tion  for  the  point  source  can  be  exploited  again  using  Eq. 
(24)  to  yield 
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An  example  of  the  fields  produced  by  a  point  source  and  a 
spiral  source  at  mid-depth  in  a  waveguide  with  h  =  20  m  is 
shown  in  Fig.  4. 

To  understand  the  differences  between  the  fields 
produced  by  the  point  source  and  the  spiral  source  in  the 
waveguide,  we  will  consider  the  expressions  for  the  fields  at 
large  distances  from  the  source,  \£np\  1.  In  this  limit, 

h£\x)  »  yTe4*-(*/ 4M™/2)]>  (53) 

and  the  expressions  for  the  point  source  and  spiral  source 
become 


In  these  expressions,  each  of  the  modes  can  be  written  as  the 
superposition  of  two  quasi-cylindrical  waves, 

and 
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propagating  at  0=±sin-1  (a Jk).  In  the  case  of  the  spiral 
source,  the  quasi-cylindrical  waves  are  reduced  in  amplitude 
relative  to  the  point  source  by 

=  cos  0,  (58) 

k 

the  same  directivity  factor  that  the  spiral  source  has  in  free 
space.  The  effect  of  this  directivity  factor  can  be  seen  in 
Fig.  4,  where  the  reverberation  in  the  upper  plot  due  to  the 
near  vertical  reflections  is  greatly  reduced  for  the  spiral 
source  in  the  lower  plot. 

The  transformation  from  the  point  source  solution  to  the 
spiral  source  solution  given  by  Eq.  (24)  can  be  applied  to 
find  the  solution  of  the  general  problem  of  a  spiral  source  in 
a  horizontally  stratified  ocean  with  an  arbitrary  sound  speed 
profile.  For  a  point  source  at  depth  z0  in  an  ocean  with  a 
sound  velocity  profile  c(z),  the  field  is  a  solution  to  the 
Helmholtz  equation, 
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where  un(x)  satisfies  the  eigenvalue  equation, 

+  -£)«.  =  o-  («i) 

To  confirm  that  the  spiral  source  solution  can  be  obtained  by 
applying  Eq.  (24)  to  Eq.  (60)  consider  the  set  of  point  sour¬ 
ces  given  by  Eq.  (1)  as  r  — >  0, 


FIG.  4.  The  transmission  loss  for  a  point  source  (upper  panel)  and  a  spiral 
source  (lower  panel),  each  located  at  z0  =  10  m  in  a  waveguide  with  a  pres¬ 
sure  release  surface  at  z  =  0  and  a  rigid  bottom  at  z  =  20  m.  In  each  example, 
/=  750  Hz  and  c  =  1500  m/s. 
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where  A0  =  —ip0ckQs.  As  the  distance  between  the  sources 
becomes  vanishingly  small,  the  singularities  can  be  written 
as 


fx  =  4nA0d-^(S(x)S(y)d(z)), 


(63) 


Similarly  for  the  pair  of  sources  along  the  y  axis  and  driven 
90°  out  of  phase, 
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In  cylindrical  coordinates, 
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and  the  sum  of  these  two  dipole  sources  becomes 
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Following  the  same  normalization  procedure  discussed  in 
Sec.  II,  the  source  term  finally  becomes 
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The  acoustic  pressure  produced  by  this  source  at  a  depth  z0 
is  described  by  the  Helmholtz  equation, 
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Assuming  a  solution  of  the  form, 
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and  integrating  both  sides  of  the  equation  with  respect  to  (j) 
from  0  to  271,  the  Helmholtz  equation  becomes 
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This  equation  can  also  be  obtained  by  taking  the  partial 
derivative  with  respect  to  p  of  both  sides  of  Eq.  (59),  which 
confirms  that  p0  is  the  point  source  solution.  The  solution  for 
a  spiral  source  in  a  waveguide  with  an  arbitrary  sound  speed 
profile  can  therefore  be  obtained  by  applying  Eq.  (24)  to 
Eq.  (60)  to  yield 
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As  an  example  of  the  application  of  the  normal  mode  solu¬ 
tions  given  by  Eqs.  (60)  and  (71),  the  pressure  fields  were 
calculated  for  sources  located  in  a  waveguide  with  a  sound 
speed  profile  similar  to  that  measured  during  the  2006  Shal¬ 
low  Water  (SW06)  experiment11  [Fig.  5(a)].  The  seafloor  in 
the  waveguide  was  a  sand  sediment  with  cb  =  1600  m/s, 
pb  =  2.0g/cm3,  and  8b  =  0.5.  The  frequency  of  each  source 
was  f—  750  Hz  and  the  normal  modes  were  calculated  using 
a  complex  mode  solver.12  The  transmission  loss  in  the  water 
column  for  each  source  out  to  3  km  is  shown  in  Figs.  5(b) 
and  5(c). 

The  fields  produced  by  the  point  and  spiral  sources  in 
Fig.  5  show  clear  differences  in  the  transmission  loss  in  the 
immediate  vicinity  of  the  sources  due  to  the  cos  8  depend¬ 
ence.  Once  the  held  has  traveled  several  water  depths  out  in 
range,  it  becomes  very  difficult  to  visually  discern  any  dif¬ 
ferences  in  the  transmission  loss.  This  can  be  seen  both  in 
Fig.  5  and  in  Fig.  6(a).  In  the  latter  plot,  the  transmission 
loss  is  plotted  at  the  source  depth  for  both  the  spiral  and  ref¬ 
erence  sources  and  these  curves  become  indistinguishable 
after  0.5  km. 


(a) 


Sound  Speed  (m/s) 


Range  (km) 


FIG.  5.  (a)  Approximation  to  the 
sound  speed  profile  measured  during 
the  2006  Shallow  Water  (SW06) 
experiment  (Ref.  11).  Using  this 
sound  speed  profile,  the  normal 
mode  solutions  were  calculated  for 
(b)  a  point  source  and  (c)  a  spiral 
source,  each  at  zo  =  30m  and  with 
/=750Hz.  For  both  solutions,  the 
bottom  was  a  sand  sediment  with 
cb=  1600 m/s,  /?6  =  2.0g/cm3,  and 
<^  =  0.5.  (d)  The  phase  difference 
between  the  two  fields  when  (/>  =  0. 
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Range  (km) 


FIG.  6.  (a)  Transmission  loss  at  the  depth  of  a  point  source  (red  line)  and  a 
spiral  source  (blue  line)  located  at  zo  =  30m  for  the  sound  speed  profile 
shown  in  Fig.  5(a).  (b)  The  phase  difference  between  a  point  source  and  a 
spiral  source  measured  at  the  source  depth. 

The  phase  difference  between  the  spiral  and  reference 
sources,  shown  in  Fig.  5(d),  is  roughly  constant  throughout 
the  water  column  for  R  >  0.5  km.  The  largest  deviations  in  the 
phase  difference  again  correspond  to  minima  in  the  propagat¬ 
ing  fields.  This  can  be  seen  more  clearly  in  Fig.  6(b),  where 
the  phase  difference  at  the  depth  of  the  sources  is  shown.  As 
the  frequency  of  the  propagating  field  is  changed,  the  loca¬ 
tions  of  the  field  minima  will  shift.  By  taking  the  mean  phase 
difference  over  the  frequency  band  of  the  transmitted  pulses, 
the  effects  of  these  nulls  can  be  mitigated  and  the  phase  dif¬ 
ference  will  be  approximately  equal  to  the  aspect  angle. 

V.  LIMITATIONS 

The  examples  and  applications  of  the  mathematical 
transformation  of  the  field  due  to  a  point  source  to  that  of  a 
spiral  source  given  by  Eq.  (24)  only  begin  to  explore  the 
tools  necessary  to  understand  the  performance  of  a  spiral 
beacon  in  an  ocean  environment.  As  more  complicated  envi¬ 
ronments  are  examined,  it  is  important  to  understand  the  lim¬ 
itations  of  this  approach.  In  Sec.  II,  solutions  to  the 
Helmholtz  equation  that  could  be  solved  as  a  sum  of  point 
sources,  such  as  using  the  method  of  images,  could  be  trans¬ 
formed  into  a  solution  for  a  collection  of  point  sources.  In 
Sec.  II,  and  again  in  Eqs.  (65)-(70)  the  point  spiral  source  is 

(a)  confined  to  the  z  axis  and  oriented  so  that  the  spiral  phase 
ramp  is  (b)  in  the  x-y  plane  and  (c)  zero  at  the  x  axis.  These 
three  conditions  are  imposed  to  achieve  the  simple  relation¬ 
ships  given  in  Eqs.  (24)  and  (26)  and  because  this  is  the 
anticipated  geometry  of  a  navigation  beacon.  For  a  collec¬ 
tion  of  sources  that  violates  (a),  one  can  simply  use  coordi¬ 
nate  transformation  and  superposition.  For  sources  that 
violate  (c),  one  needs  to  only  change  the  phase  term  in  Eqs. 
(24)  and  (71)  and  again  invoke  superposition.  Violation  of 

(b)  is  more  complicated,  but  can  also  be  accomplished  with 
careful  coordinate  transformation  and  superposition. 

Consider  the  simple  case  of  a  wedge  environment  with  a 
sloping,  pressure  release  bottom  and  a  flat  pressure  release 
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surface.13  Again  this  problem  can  be  solved  with  the  method 
of  images,  however  these  images  no  longer  lie  on  the  z  axis. 
In  order  to  satisfy  the  boundary  conditions  for  the  surface, 
the  image  source  will  again  be  placed  on  the  z  axis  but  for 
the  sloping  bottom,  the  image  source  must  lie  along  a  line 
that  is  normal  to  the  bottom  and  hence  no  longer  parallel  to 
the  z  axis.  The  spiral  source  solution  to  this  environment  can 
be  constructed  from  image  sources  as  well,  but  the  plane  of 
these  sources  must  the  tilted  in  order  to  match  the  boundary 
conditions  and  hence  Eq.  (24)  cannot  be  applied  directly. 

This  requirement  that  the  environment  be  horizontally 
stratified  and  range  independent  becomes  important  as  well 
when  the  effects  of  roughness  or  volume  scattering  are  con¬ 
sidered.  These  effects  are  typically  addressed  for  propaga¬ 
tion  from  a  point  source  through  approximate  methods,  such 
as  perturbation  theory,14  or,  more  recently,  using  numerical 
techniques  such  as  rough  surface  parabolic  equation  codes15 
or  finite  difference  time-domain  methods.16  While  the  theo¬ 
retical  approximations  typically  find  plane  wave  solutions 
and  hence  can  be  applied  to  the  spiral  source  through  the 
techniques  given  in  Sec.  Ill,  the  numerical  approaches  solve 
a  form  of  the  Helmholtz  equation  at  the  boundary  or  in  the 
volume  for  the  full  field  from  the  point  source.  Since  these 
rough  surfaces  involve  facets  of  the  boundary  that  are  no 
longer  horizontal,  the  transformation  given  by  Eq.  (24)  can¬ 
not  be  applied  to  the  scattered  field.  In  order  to  apply  numer¬ 
ical  techniques  such  as  finite  element  solutions,  the  full 
“spiral”  Helmholtz  equation  given  by  Eq.  (68)  would  need  to 
be  solved. 

VI.  CONCLUSIONS 

In  order  to  understand  the  performance  of  a  spiral  source 
and  spiral  beacon  in  an  ocean  environment,  the  spiral  analog 
of  the  acoustic  point  source  was  developed.  This  ideal  spiral 
source  was  constructed  from  two  perpendicular  dipoles 
placed  at  the  origin  and  driven  90°  out  of  phase.  The  result¬ 
ing  source  produces  a  pressure  field  that  has  the  required 
phase  dependence  for  a  spiral  wave  front.  This  spiral  source 
was  shown  to  be  related  through  a  simple  transform  to  the 
point  source.  The  transformation  was  applied  to  determine 
the  field  propagating  from  a  spiral  source  near  an  ocean  sur¬ 
face  or  bottom  and  in  a  horizontally  stratified  waveguide. 
This  solution  technique  is  limited  to  range-independent 
problems,  but  the  spiral  wave  front  Helmholtz  equation 
given  in  Sec.  IV  could  be  used  to  develop  numerical 
approaches,  such  as  parabolic  equation  codes,  to  solve  for 
the  spiral  field  directly  in  range-dependent  environments. 
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Transducers  for  acoustic  beacons  which  can  produce  outgoing  signals  with  wave  fronts  whose  hori¬ 
zontal  cross  sections  are  circular  or  spiral  are  studied  experimentally.  A  remote  hydrophone  is  used 
to  determine  its  aspect  relative  to  the  transducers  by  comparing  the  phase  of  the  circular  signal  to 
the  phase  of  the  spiral  signal.  The  transducers  for  a  “physical-spiral”  beacon  are  made  by  forming  a 
strip  of  1-3  piezocomposite  transducer  material  around  either  a  circular  or  spiral  backing.  A 
“phased- spiral”  beacon  is  made  from  an  array  of  transducer  elements  which  can  be  driven  either  in 
phase  or  staggered  out  of  phase  so  as  to  produce  signals  with  either  a  circular  or  spiral  wave  front. 
Measurements  are  made  to  study  outgoing  signals  and  their  usefulness  in  determining  aspect  angle. 
Vertical  beam  width  is  also  examined  and  phase  corrections  applied  when  the  hydrophone  is  out  of 
the  horizontal  plane  of  the  beacon.  While  numerical  simulations  indicate  that  the  discontinuity  in 
the  physical- spiral  beacon  introduces  errors  into  the  measured  phase,  damping  observed  at  the  ends 
of  the  piezocomposite  material  is  a  more  significant  source  of  error.  This  damping  is  also  reflected 
in  laser  Doppler  vibrometer  measurements  of  the  transducer’s  surface  velocity. 
[http://dx.doi.Org/10.1121/l.3699170] 
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I.  INTRODUCTION 

A  technique  for  underwater  navigation  has  recently  been 
developed  which  utilizes  a  beacon  capable  of  producing  a  spi¬ 
ral  wave  front.1,2  The  beacon  consists  of  two  cylindrical  sour¬ 
ces,  lying  in  the  x-y  plane,  that  radiate  sound  from  their 
circumferential  surfaces.  One,  a  “reference  source”  produces 
a  wave  front  with  a  cross  section  in  the  x-y  plane  that  forms 
concentric  circles.  A  second  transducer,  a  “spiral  source,”  pro¬ 
duces  a  wave  front  with  a  cross  section  that  forms  a  linear  spi¬ 
ral.  A  remote  hydrophone  can  then  determine  its  aspect 
relative  to  the  beacon  by  comparing  the  phase  of  the  two  sig¬ 
nals.  Theoretical  groundwork  for  the  spiral  wave  front  beacon 
is  presented  in  a  companion  paper2  and  the  reader  is  referred 
there  for  a  complete  discussion  of  the  concept. 

The  concept  of  a  spiral  wave  field  has  been  previously 
investigated  by  Ceperley.3  A  spiral  wave  front  can  be  pro¬ 
duced  by  either  a  single  element  transducer  in  the  shape  of  a 
linear  spiral,  a  “physically  induced”  spiral  wave  front,  or  by 
using  an  array  of  phased  elements,  a  “phased”  spiral  wave 
front.  As  previously  shown,2  the  properties  of  each  of  these 
wave  fronts  is  slightly  different  and  each  has  advantages  and 
disadvantages.  This  paper  discusses  the  design,  construction, 
and  testing  of  prototypes  of  both  a  physically  induced  and  a 
phased  spiral  wave  front  beacon. 

Details  of  the  prototype  reference  and  spiral  sources  are 
given  in  Sec.  II.  Each  source  is  tested  in  both  underwater 
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laboratory  facilities  and  in  air  using  laser  Doppler  vibrome- 
try.  Details  of  the  testing  are  given  in  Sec.  Ill  followed  by  a 
brief  theoretical  discussion  in  Sec.  IV  of  how  aspect  is  deter¬ 
mined  by  comparison  of  the  reference  and  spiral  sources. 
This  discussion  is  required  to  interpret  the  in-plane  and  out- 
of-plane  results  that  follow  in  Sec.  V.  It  is  shown  that  that 
spiral  wave  fields  can  be  realized  by  using  either  of  the 
physically  induced  or  phased  techniques. 

II.  TRANSDUCER  DESIGN  AND  CONSTRUCTION 

Figure  1  is  a  photograph  of  the  physical  beacon  with  the 
physically  induced  spiral  wave  front  transducer.  The  refer¬ 
ence  source  is  stacked  above  the  spiral  source.  On  the  spiral 
source,  the  one  wavelength  spiral  offset  is  located  near  the 
bottom  of  the  photograph.  The  beacon  was  manufactured  by 
Dr.  Thomas  Howarth  and  Kim  Benjamin  at  NAVSEA  New¬ 
port  Division  in  Newport,  RI.  To  construct  the  transducers,  a 
6.35  mm  thick  strip  of  1-3  piezocomposite  (Navy  type  I) 
transducer  material  was  formed  around  a  sound  absorbing 
backing  material.  After  attaching  leads  to  the  transducer,  a 
layer  of  polyurethane  was  formed  in  a  thin  layer  over  the 
piezocomposite  and  the  backing  as  a  water  encapsulant.  The 
height  of  each  active  surface  is  19.1  mm  and  the  overall 
height,  including  the  outer  polyurethane  coating,  is  29.8  mm. 
The  reference  source  has  an  outer  radius  of  82.6  mm.  This  is 
also  the  average  radius  of  the  spiral  source,  a^.  However,  the 
radius  for  the  spiral  source  varies  linearly  from  73.0  mm  to 
93.1  mm  giving  an  offset  of  19.1  mm,  equal  to  the  opera¬ 
tional  wavelength.  In  addition,  fixtures  were  added  so  that 
the  transducers  can  be  arranged  to  form  a  beacon  as  seen  in 
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FIG.  1.  Photograph  of  spiral  wave  front  beacon  prototype.  The  circular  ref¬ 
erence  source  is  at  the  top  and  the  spiral  source  is  at  the  bottom.  The  discon¬ 
tinuity  in  the  transducer  material  and  the  spiral  offset  are  indicated. 


Fig.  1.  A  hollow  steel  rod  protruding  from  the  top  of  the  bea¬ 
con  supports  the  beacon  and  keeps  it  rigidly  fixed  to  the  rota¬ 
tional  stage  during  the  measurements.  This  rod  also  houses 
the  wires  powering  the  transducers.  When  the  beacon  is  sus¬ 
pended  for  testing,  the  spiral  source  is  below  the  x-y  plane  at 
—A z,  and  the  distance  between  the  vertical  centers  of  the 
transducers  is  2 A z  =  31.9  mm.  Each  of  the  transducers  has  a 
discontinuity  where  the  ends  of  the  composite  material  come 
together.  For  the  spiral  source,  this  is  located  at  the  offset 
which  is  defined  as  0°  aspect  angle.  The  discontinuity  in  the 
reference  transducer  is  located  at  90°.  Both  the  discontinuity 
and  the  offset  are  indicated  in  Fig.  1.  The  configuration  of 
the  beacon  and  the  coordinate  system  are  shown  in  Fig.  2. 

A  phased  transducer  array  was  also  constructed  to  be 
used  as  a  spiral  source.  The  array  is  similar  in  design  to  the 
circular  reference  source  in  Fig.  1  and  has  the  same  dimen¬ 
sions.  However,  the  electrical  contact  material  is  etched  so 
that  portions  of  the  3-1  piezocomposite  transducer  material 
can  act  as  separate  elements.  There  are  16  equally  spaced 


Y 


FIG.  2.  Beacon  configuration  and  orientation  showing  aspect  angle,  </>,  and 
tilt  angle,  6,  to  an  observer  at  point  P.  Also  indicated  is  the  vertical  distance 
between  the  centers  of  the  transducers,  2 A z,  and  the  average  radius  of  the 
transducers,  ao,  measured  from  the  z-axis.  The  support  rod  is  shown  protrud¬ 
ing  from  the  top  of  the  circular  transducer. 


elements  each  covering  22.5°.  Due  to  manufacturing  error, 
one  pair  of  neighboring  elements  are  shorted  and  act  as  a  sin¬ 
gle  45°  element. 

The  operational  frequency  of  the  physically  induced  ver¬ 
sion  of  the  beacon  is  such  that  the  offset  distance,  19.1  mm, 
corresponds  to  one  wavelength.  The  water  temperature  in  the 
acoustic  test  facility  at  NAVSEA  Panama  City  Division  in 
Panama  City,  FL,  where  the  experiments  are  carried  out,  var¬ 
ied  between  28.0 °C  and  28.1  °C  on  the  test  dates.  The  water 
in  the  test  facility  is  fresh  and  the  sound  speed  is  calculated4 
to  be  1503  m/s  corresponding  to  an  operational  frequency  of 
78.9  kHz.  Ideally,  the  sound  speed  would  be  measured  using 
time  of  flight  measurements.  However,  the  error  introduced 
by  incorrectly  calculating  the  sound  speed  is  likely  small  as 
an  error  of  10  m/s  would  result  in  an  aspect  error  of  only  2.4° 
at  the  offset.  Experiments  involving  the  phased  transducer 
array  were  carried  out  in  a  freshwater  tank  at  the  Naval 
Research  Faboratory  in  Washington,  DC.  The  sound  speed 
is  not  a  factor  in  the  operation  of  the  phased  spiral  transducer 
array. 

III.  MEASUREMENTS  AND  SIGNAL  PROCESSING 
A.  Acoustic  testing 

The  measurements  of  each  transducer  were  performed 
using  a  simple  send  and  receive  configuration.  The  beacons 
are  suspended  in  the  water  tank  by  a  computer  controlled 
rotational  stage.  Each  transducer  on  the  beacon  was  driven 
independently  through  its  own  power  amplifier.  The  hydro¬ 
phone  was  positioned  1.2  m,  approximately  62  wavelengths, 
from  the  center  of  the  beacon.  Aspect  angle,  </>,  was  varied  in 
increments  of  1°.  The  tilt  angle,  0,  was  varied  by  raising  or 
lowering  the  hydrophone  out  of  the  plane  of  the  beacon.  The 
results  were  adjusted  for  spherical  spreading  to  1  m  from  the 
center  of  the  beacon  to  the  hydrophone  position,  which  was 
determined  by  time-of-flight. 

The  transducers  in  the  physical  spiral  beacon  are  driven 
with  a  series  of  three  tone  bursts.  First,  a  tone  burst  was  sent 
from  the  reference  transducer,  then  from  the  spiral  transducer, 
and  then  from  both  together.  Each  tone  burst  is  20  cycles  with 
a  cosine  squared  window  over  the  first  and  last  10%  of  the  sig¬ 
nal.  A  repetition  rate  of  10  Hz  is  sufficient  for  the  dissipation 
of  reverberations  in  either  test  facility  at  these  frequencies. 
The  signals  were  received  on  a  Navy  standard  H-52  hydro¬ 
phone5  whose  free-field  voltage  response  at  78.9  kHz  is  cali¬ 
brated  at  —178.4  dB  re  1  V//iPa  at  the  hydrophone  face.  The 
received  signal  was  first  band-pass  filtered  between  39  kHz 
and  158  kHz.  Each  of  the  three  tone  bursts  were  then  isolated 
by  identifying  the  maximum  of  a  cross-correlation  and  win¬ 
dowing  away  the  rest  of  the  signal.  Since  these  are  not  pure 
continuous  wave  (CW)  signals,  care  must  be  taken  when 
determining  the  phase  of  the  tone  bursts.  The  Fourier  trans¬ 
form  of  each  tone  burst  was  interpolated  to  78.9  kHz  to  calcu¬ 
late  the  amplitude  A  and  the  phase  ®.  Each  of  the  three  tone 
bursts  is  written  as  a  CW  signal  with  a  complex  amplitude, 
Acel®c  for  the  reference  return,  Asel®s  for  the  spiral  return,  and 
At,elQ>b  for  the  return  when  both  transducers  are  activated  to¬ 
gether.  This  puts  the  data  into  a  convenient  format  to  measure 
the  phase  difference  and  to  calculate  the  aspect. 
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Testing  with  the  phased  array  was  performed  using  only 
a  single  windowed  tone  burst  without  a  corresponding  refer¬ 
ence  transducer  since  the  array  can  be  used  for  both  the  spiral 
and  reference  signals.  Because  the  phased  array  does  not  rely 
on  a  physical  gap  to  produce  a  spiral  wave  front  it  was  tested 
with  26.3  kHz,  52.6  kHz,  79.0  kHz,  and  105.3  kHz  tone  bursts 
of  10,  20,  30,  and  40  cycles,  respectively.  A  Reson  TC  4013 
hydrophone  (Slangerup,  Denmark)  with  a  receiving  sensitivity 
of  —211  dB  re  1  V/^Pa  at  1  m  was  used  for  data  collection. 
With  the  hydrophone  in  plane,  0  =  0°,  the  output  from  each 
element  was  measured  independently.  Otherwise,  the  meas¬ 
urements  were  conducted  identically  to  that  of  the  physical 
beacon’s  transducers.  The  outputs  of  the  phased  transducer 
array  as  a  reference  and  spiral  source  were  then  synthesized 
from  the  results. 

B.  Laser  Doppler  vibrometer  testing 


hydrophone  is  moved  out-of-plane,  0^0,  the  difference 
between  the  phase  of  the  spiral  transducer  and  the  reference 
varies  as  a  function  of  their  separation  distance.  For  the  anal¬ 
ysis  here,  however,  a  smooth  phase  ramp  will  be  assumed 
for  in-plane  measurements  with  appropriate  corrections 
made  for  out-of-plane  measurements.  Deviations  from  ideal 
behavior  due  to  diffractive  effects  will  not  be  taken  into 
account,  but  comparisons  between  the  experimental  and  nu¬ 
merical  results  are  examined. 

With  the  hydrophone  in-plane,  0^0,  the  incoming  sig¬ 
nal  from  the  reference  transducer  has  phase 

O  c  =  k(R-ao),  (1) 

where  R  is  the  distance  to  the  center  of  the  transducer  and  ao 
is  the  radius  of  the  transducer.  The  phase  of  the  signal  from 
the  spiral  transducer  has  a  more  complicated  form  due  to  the 
physical  phase  ramp.  It  can  be  written 


In  addition  to  acoustic  tests,  the  surface  velocities  of  the 
physical  spiral  transducers  in  air  were  also  measured  using  a 
laser  Doppler  vibrometer6  (LDV)  at  the  Naval  Research  Labo¬ 
ratory  in  Washington,  DC.  The  surface  of  the  transducer  was 
painted  with  a  highly  reflecting  white  paint  and  mounted  on  a 
rotational  stage  sitting  on  a  concrete  floor  for  stability.  The 
laser  is  focused  onto  a  spot  in  the  center  of  the  transducer  and 
the  transducer  rotated  between  —180°  and  180°  at  increments 
of  1°.  The  transducer  is  driven  with  a  100  cycle  tone  burst  and 
the  output  from  the  power  amplifier  and  the  LDV  recorded 
simultaneously  so  their  relative  phase  could  be  measured.  The 
LDV  also  generated  an  additional  direct  current  voltage  pro¬ 
portional  to  the  strength  of  the  reflected  signal.  This  signal, 
recorded  simultaneously,  is  used  to  reject  data  points  which 
have  inadequate  signal  levels. 

Twenty  tone  bursts  at  each  position  were  recorded  and 
averaged.  The  phase  and  amplitude  of  the  signals  were  cal¬ 
culated  from  the  Fourier  transform  of  the  signals  at  the  driv¬ 
ing  frequency  78.9  kHz.  Because  the  surface  of  the 
transducers  is  not  smooth,  the  results  are  quite  noisy.  Also, 
since  the  active  material  is  a  1-3  piezocomposite,  there  are 
rods  of  active  piezoceramic  imbedded  in  a  lossy  epoxy  sub¬ 
strate.  As  the  laser  moves  over  the  surface  it  is  also  moving 
over  these  rods,  leading  to  variation  in  amplitude.  To  mini¬ 
mize  these  effects,  a  sliding  average  is  applied  to  the  data 
where  each  data  point  is  replaced  by  the  average  of  itself  and 
two  data  points  to  either  side.  For  the  case  where  a  data  point 
is  dropped  due  to  low  signal  levels,  the  remaining  samples 
are  averaged.  The  results  are  then  converted  to  velocities 
using  the  calibration,  20  mm/(s  V),  provided  by  the  vibrome¬ 
ter’ s  manufacturer. 


$ s  =  k{R-ao)-(i)  +  n A  (2) 

m 

for  —  n  <  (j)  <  7i.  The  correction,  </>/|</>|,  gives  the  sign  of  </> 
and  accounts  for  the  discontinuity  at  0  =  0.  Note  that  since 
the  coordinate  </>  is  circular,  the  equations  could  have  been 
written  with  </>  running  between  0  and  2n  without  the  phase 
correction.  However,  this  convention  puts  the  spiral  offset  at 
the  left  and  right  edges  of  subsequent  figures  and  makes 
analysis  more  difficult.  The  aspect,  </>,  to  the  target  can  then 
be  calculated  by  simply  subtracting  the  phases  of  the 
returned  signals 

0  =  A<5  +  7tA,  (3) 
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where  AO  =  Oc  —  O^.  Using  symmetry,  Eq.  (3)  can  be 
rewritten  as 


</> 


=  AO  —  n 


AO 

W 


(4) 


If  only  the  amplitudes  of  the  signals  are  available,  the 
aspect  can  be  calculated  using  the  combined  return  as  well. 
The  superposition  of  the  spiral  and  reference  returns  can  be 
written  as 


Ahe‘^  =  Acel(^  +  (5) 

Taking  the  square  of  the  absolute  value  of  Eq.  (5)  and  insert¬ 
ing  Eqs.  (1)  and  (2)  for  the  phase,  the  aspect  can  be  deter¬ 
mined  from 


IV.  DETERMINATION  OF  ASPECT 

Design  of  a  spiral  transducer  ideally  gives  a  smooth 
phase  ramp  with  aspect,  </>,  at  the  operational  frequency. 
However,  the  wave  field  generated  by  the  real  transducer  is 
complex  and  the  phase  ramp  is  not  perfectly  smooth  due  to 
the  spiral  offset  and  proximity  to  the  reference  transducer,2 
in  the  case  of  the  physical  transducer,  or  the  finite  aperture 
and  side  lobe  interference  in  the  case  of  the  phased  trans¬ 
ducer  array.  In  addition,  with  the  physical  beacon,  as  the 


(j)  =  cos  1 


A2  —  A2  —  A2 


2 ArAs 


</> 

+*m 


1.6) 


While  Eq.  (6)  requires  only  the  amplitudes,  apart  from  the 
phase  correction,  it  is  ambiguous  since  the  inverse  cosine  is 
multivalued.  For  this  reason,  Eq.  (4)  will  be  used  to  calculate 
the  results  presented  here.  Note  that  when  using  either  of 
these  methods,  the  results  do  not  need  to  be  “unwrapped”  as 
the  majority  of  the  phase,  k(R  +  ao),  is  subtracted  out. 


3750  J.  Acoust.  Soc.  Am.,  Vol.  1 31 ,  No.  5,  May  201 2 


B.  R.  Dzikowicz  and  B.  T.  Hefner:  Spiral  wave  front  beacon 


Downloaded  16  Jul  2012  to  128.95.76.64.  Redistribution  subject  to  ASA  license  or  copyright;  see  http://asadl.org/journals/doc/ASALIB-home/info/terms.jsp 


When  the  hydrophone  is  out-of-plane  for  the  physical 
beacon,  a  far-held  correction  can  be  applied  as  given  by  Hef¬ 
ner  and  Dzikowicz,2 


</> 


AO 
cos  8 


471  Az 


tan  8  —  n 


AO 

m 


where  Az  is  as  described  in  Sec.  II  and  Fig.  2. 


(7) 


V.  RESULTS  OF  THE  TESTING 

A.  In-plane  results  for  the  physically  induced  spiral 

The  aspect  angle  determined  from  the  physically 
induced  beacon  at  8  =  0°,  is  shown  in  Fig.  3.  The  result  is 
calculated  using  the  data  collection  and  signal  processing 
discussed  in  Sec.  Ill  and  Eq.  (4).  For  the  most  part,  the 
results  follow  the  position  of  the  rotational  stage  well,  except 
at  4>  7^  0°,  the  position  of  the  offset  in  the  spiral  transducer, 
and  at  </>  =  90°,  the  position  of  the  discontinuity  in  the  refer¬ 
ence  transducer  described  in  Sec.  II  and  shown  in  Fig.  1.  The 
remainder  of  the  results  and  discussion  will  be  presented 
with  these  deviations  in  mind  and  an  eye  toward  understand¬ 
ing  their  cause. 

The  amplitudes  of  the  tone  bursts  for  the  in-plane, 
8  =  0°,  case  are  shown  in  Fig.  4.  Although  there  is  large  var¬ 
iation  in  the  response  with  aspect,  the  expected  behavior  is 
observed:  a  flat  response  for  either  of  the  transducers  fired 
individually,  but  an  amplitude  which  varies  with  respect  to 
aspect,  when  both  are  fired  together.  Dips  in  the  signal 
appear  both  in  the  reference  transducer  amplitude  near  the 
discontinuity  and  in  the  spiral  transducer  at  the  position  of 
the  spiral  offset.  The  source  levels  of  each  received  tone 
burst  averaged  over  all  aspects  are  165,  166,  and  168  dB  re 
1  ,uPa  at  1  m.  Note  that  apart  from  the  position  of  the  discon- 


FIG.  3.  Position  as  determined  from  the  physical  beacon  data  using  Eq.  (4) 
compared  to  the  true  position  of  the  rotational  stage.  The  largest  deviations 
from  true  are  at  (j>  =  0°,  the  position  of  the  spiral  offset,  and  at  (j)  =  90°, 
where  there  is  a  discontinuity  in  the  reference  transducer.  These  positions 
are  indicated  with  vertical  gray  lines. 


180° 


180° 


FIG.  4.  Amplitude  of  each  tone  burst  with  respect  to  rotator  position.  The 
top  plot  corresponds  to  the  first  hydrophone  return  as  shown  in  Fig.  4.  The 
middle  and  bottom  plots  correspond  to  the  second  and  third  hydrophone 
returns  in  Fig.  4,  respectively.  The  vertical  line  at  <f>  =  0°  corresponds  to  the 
spiral  offset  and  the  vertical  line  at  (j>  =  90°  corresponds  to  the  position  of 
discontinuity  in  the  reference  transducer,  as  shown  in  Fig.  1. 


tinuity  and  the  offset,  the  calculated  aspect  angle  in  Fig.  3 
does  not  vary  as  the  amplitudes  shown  in  Fig.  4.  In  addition, 
by  taking  the  phase  difference  there  is  no  correction  needed 
for  the  range  to  the  beacon.  For  example,  a  wobble  in  the 
support  rod  between  the  rotational  stage  and  the  beacon  is 
revealed  as  a  large,  approximately  cosine  oscillation  in  the 
phase  of  both  transducers.  The  wobble  is  about  10  mm  off 
axis  at  the  beacon  and  is  likely  a  result  of  the  rod  being  bent 
or  not  attached  perfectly  normal  to  the  rotational  stage.  This 
phase  shift  due  to  the  wobble  cancels  out  in  the  calculation 
and  is  not  seen  in  the  results  in  Fig.  3. 

Numerical  results  using  the  finite  element  method  (FEM), 
shown  in  Fig.  9  of  Ref.  2,  predict  variation  in  the  phase  due  to 
the  geometry  of  the  spiral  offset.  The  numerical  results  repre¬ 
sent  an  absolute  limit  to  the  accuracy  of  the  prototype  beacon. 
For  practical  operation,  it  may  be  possible  to  calibrate  the 
transducer  and  eliminate  these  effects.  Figure  5  reproduces 
these  results  and  compares  them  with  the  difference  between 
the  calculated  results  and  the  ideal  results  in  Fig.  3.  There  is 
little,  if  any,  agreement  between  the  numerical  results  and  the 
experimental.  Other  effects  at  the  offset  in  the  spiral  trans¬ 
ducer  and  the  discontinuity  in  the  reference  transducer  are 
likely  overshadowing  any  effect  from  the  geometry  itself. 

In  order  to  understand  the  dips  in  the  amplitude  and  the 
error  in  the  angle  calculation  at  the  spiral  offset  in  the  spiral 
transducer  and  at  the  discontinuity  in  the  reference  transducer, 
the  amplitudes  are  compared  to  the  surface  vibration  ampli¬ 
tudes  as  measured  by  the  LDV  as  described  in  Sec.  Ill  B. 
These  comparisons  are  shown  in  Fig.  6.  The  acoustic  drop  in 
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Rotator  Position,  <j> 

FIG.  5.  Phase  difference  between  the  calculated  phase  and  the  rotator 
position  in  Fig.  5  compared  with  the  numerical  (FEM)  calculation  of  the 
phase  error  based  on  the  geometry  of  the  prototype.  The  vertical  line  at 
<p  =  0°  corresponds  to  the  spiral  offset  and  the  vertical  line  at  (p  =  90°  cor¬ 
responds  to  the  position  of  discontinuity  in  the  reference  transducer,  as 
shown  in  Fig.  1. 

the  signal  levels  at  the  offset  and  discontinuity  is  verified  by 
similar  drops  in  the  LDV  data.  Note  that  the  sharp  dip  that 
appears  on  all  three  tone  bursts  at  0  ^  —155°  in  both  Figs.  4 
and  6  is  not  reflected  in  the  LDV  data,  nor  does  it  appear  in 
the  unaveraged  LDV  data.  This  indicates  that  this  is  more 
likely  a  drop-out  or  glitch  in  the  electronics  at  that  position 
during  the  acoustic  data  collection. 

A  more  revealing  result  is  the  comparison  shown  in 
Fig.  5  of  the  phase  differences  in  the  LDV  data  to  the  phase 
error  in  the  acoustic  data.  Many  of  the  features  match  up, 
especially  at  the  offset  and  discontinuity.  This  indicates  that 
errors  in  phase  come  directly  from  the  vibration  of  the  trans¬ 
ducer  itself  rather  than  from  experimental  or  signal  process¬ 
ing  problems.  Since  the  LDV  data  is  collected  in  air  and  the 
acoustic  data  in  water,  this  result  also  implies  that  the  load¬ 
ing  from  the  water,  in  addition  to  that  of  the  polyurethane 
layer  present  in  both  experiments,  does  not  strongly  affect 
the  phase  difference  between  the  transducers,  another  indica¬ 
tion  of  the  robustness  of  this  technique. 

B.  Out  of  plane  results  for  the  physically  induced 
spiral 

This  prototype  beacon  is  designed  to  have  a  narrow 
beam  pattern.  The  amplitude  as  a  function  of  angle  for  five 
different  pitch  angles  is  shown  in  Fig.  7.  The  experimental 
beam  width  suggested  by  these  measurements  is  narrower 
than  is  expected  from  finite  element  analysis.  Recall  that  the 
reference  transducer  is  above  the  spiral  transducer  and  when 
180°  >  (j)  >  0°,  the  radius  of  spiral  transducer  is  smaller 
than  that  of  the  reference  transducer.  Thus,  when  9  >  0°,  the 
signal  from  the  spiral  transducer  can  be  partially  obscured 
by  the  reference  transducer.  The  opposite  is  true  when  below 
the  plane  or  on  the  opposite  side  of  the  beacon.  However,  if 


Rotator  Position,  <p 

FIG.  6.  Acoustic  amplitude  of  the  tone  bursts  from  each  of  the  transducers 
compared  with  the  surface  velocity  as  determined  using  laser  Doppler  vibr- 
ometry  in  air.  The  vertical  line  at  <p  =  0°  corresponds  to  the  spiral  offset  and 
the  vertical  line  at  (p  =  90°  corresponds  to  the  position  of  the  discontinuity 
in  the  reference  transducer,  as  shown  in  Fig.  1. 

this  has  any  effect  at  the  angles  tested,  it  is  obscured  by  other 
effects  when  looking  at  the  unaveraged  amplitude  of  the 
signal. 

In  order  to  calculate  aspect  angle  when  9  ^  0°,  the  cor¬ 
rection  given  in  Eq.  (7)  must  be  applied.  Figure  8  shows  the 
effect  of  this  correction  for  two  of  the  pitches  in  Fig.  6,  —25° 


FIG.  7.  Mean  peak  pressure  of  each  transducer  as  a  function  of  pitch  angle 
suggesting  an  approximate  vertical  beamwidth. 
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FIG.  8.  Application  of  the  correction  in  Eq.  (7)  to  out-of-plane  data.  The 
top  line  represents  the  uncorrected  data  at  6  =  —25°,  and  the  bottom  line, 
uncorrected  data  at  6  =  23°.  Arrows  indicate  the  effect  of  the  correction. 
The  vertical  line  at  =  0°  corresponds  to  the  spiral  offset  and  the  vertical 
line  at  <f>  =  90°  corresponds  to  the  position  of  discontinuity  in  the  reference 
transducer,  as  shown  in  Fig.  1. 

and  23°.  Prior  finite  element  work  indicates  that  when  the 
distance  between  the  transducers  is  small,  their  signals  inter¬ 
fere  with  each  other.2  This  is  especially  true  for  measure¬ 
ments  out-of-plane.  However,  despite  this  interference  and 
the  lower  signal  level,  the  results  after  correction  are  reason¬ 
able  and  validate  the  correction  given  in  Eq.  (7). 

C.  Phased  array  results 

For  the  phased  array  the  output  from  each  element  is 
measured  separately.  The  response  of  a  single  element  is 
shown  in  Fig.  9.  Two  lines  denoted  with  “A”s  show  the 
physical  angular  extent  of  the  element  on  the  array.  Note 
that  at  these  frequencies,  the  spreading  of  the  beam  is  pri¬ 
marily  a  result  of  the  curvature  of  the  element  itself  rather 


FIG.  9.  Polar  plot  of  the  pressure  amplitude  from  a  single  element  in  the 
phased  transducer  at  several  frequencies.  The  lines  denoted  by  “A”  show  the 
physical  angular  extent  of  the  element. 
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FIG.  10.  Polar  plot  of  the  pressure  amplitude  from  each  element  at  78.9 
kHz.  Elements  8  and  9  (counting  counterclockwise  from  the  x-axis)  are 
shorted  together  and  thus  act  as  a  wider  element.  The  center  of  this  wider 
element  is  marked  with  a  gray  line  corresponding  to  the  angle  marked  by  a 
similar  gray  line  in  Fig.  11. 

than  its  aperture.  Figure  10  shows  the  response  from  each 
element  at  78.9  kHz.  A  manufacturing  error  caused  elements 
8  and  9  (counting  counterclockwise  from  the  x-axis)  to  be 
shorted  together  within  the  polyurethane  coating.  Element 
capacitance  and  FDV  measurements,  too  mundane  to  other¬ 
wise  report  here,  support  this.  Since  activation  of  either 
element’s  lead  drives  both  elements,  there  is  a  large  contri¬ 
bution  at  this  position  further  denoted  by  a  gray  line  above 
the  —x-axis. 

In  order  to  test  the  accuracy  of  the  phased  beacon  and 
compare  the  results  to  the  physical  beacon,  the  reference  and 
spiral  signals  are  synthesized  by  the  superposition  of  the  sin¬ 
gle  element  results  at  78.9  kHz.  For  the  reference  signal,  the 
signals  from  each  element  are  added  in  phase.  For  the  spiral 
signal,  a  phase  shift  is  added  to  each  element  corresponding 
to  its  position  on  the  array.  In  both  cases,  the  contribution  of 
the  shorted  elements  is  halved  since  each  element  is  being 
driven  twice.  Aspect  determination  is  simulated  by  subtract¬ 
ing  the  phase  of  these  two  signals  as  in  Eq.  (4).  The  results 
are  shown  in  Fig.  11.  During  the  course  of  calculation,  a  simi¬ 
lar  wobble  in  the  support  rod,  as  found  with  the  physical  bea¬ 
con  and  discussed  in  Sec.  V  A,  is  discovered.  However,  as 
before,  when  the  reference  signal  is  subtracted,  this  vanishes. 

VI.  DISCUSSION 

Figures  3  and  1 1  show  that  the  spiral  beacon  concept  is 
certainly  feasible.  Comparison  of  calculated  aspect  to  true 
aspect  gives  a  root  mean  squared  error  of  9.1°  for  the  physi¬ 
cal  beacon  and  10.1°  for  the  phased  beacon. 

Results  from  the  FDV  experiments  on  the  physical  bea¬ 
con  show  that  the  biggest  errors,  at  the  spiral  offset  and  at 
the  reference  dislocation,  are  not  due  to  the  geometry  of  the 
system,  but  rather,  are  related  to  discontinuities  where  the 
ends  of  the  strips  of  active  material  come  together.  There  is 
a  relationship  between  the  reduced  amplitude  and  the  phase 
shifts  at  these  positions.  It  can  be  seen  by  examining 
Figs.  3-5  that,  where  the  strips  come  together,  there  is  a 
decrease  in  both  amplitude  and  phase  for  both  the  reference 
and  spiral  transducers.  Treating  the  surface  of  the  transducer 
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FIG.  1 1 .  Position  as  determined  from  the  acoustic  beacon  data  using  Eq.  (4) 
compared  to  the  true  position  of  the  rotational  stage.  f  =  —11.25°,  is  the 
position  directly  between  the  two  shorted  elements.  This  position  is  indi¬ 
cated  with  a  vertical  gray  line. 


as  a  damped-driven  harmonic  oscillator  allows  for  simple 
qualitative  analysis.7  Defining  the  mechanical  impedance  as 
Z  =  R  +  iX  where  R  and  X  are  both  real,  the  velocity  and 
phase  of  the  oscillator  are 


\u\ 


F 

M 


and 


O  =  tan 


X 

R 


(8) 


different  transducers.  Future  designs  may  align  the  disconti¬ 
nuities,  minimizing  the  effect  of  damping. 

The  overall  error  for  the  phased  beacon  is  comparable  to 
that  of  the  physical  beacon.  Although  the  physical  beacon 
requires  far  simpler  electronics  to  implement,  there  are  two 
advantages  for  using  the  phased  beacon  approach.  First,  the 
device  can  operate  at  a  range  of  frequencies  or  even  at  multi¬ 
ple  frequencies.  Second,  there  is  no  out-of-plane  or  vertical 
correction  to  make  since  the  reference  signal  and  the  spiral 
signal  are  sent  from  the  same  set  of  elements. 

Lessons  from  these  tests  show  that  even  small  levels  of 
uneven  mechanical  damping  can  cause  unwanted  variation 
in  the  phase  and  thus  the  accuracy  of  the  device.  Future 
designs  must  consider  these  effects.  Damping  considerations 
aside,  designing  a  commercial  beacon  then  becomes  a  matter 
of  making  trade-offs  on  several  parameters.  Narrow  vertical 
beam  widths  might  be  best  for  long  ranges  or  when  the 
hydrophone  and  beacon  are  at  known  depths,  and  wide  verti¬ 
cal  beamwidths  are  preferred  in  close  range  applications 
such  as  docking.  Large  vertical  spacing  between  the  trans¬ 
ducers  provides  a  smoother  phase  shift  in  aspect,  but  requires 
greater  correction  when  pitch  varies.  Also,  the  frequency 
must  be  chosen  so  as  not  to  interfere  with  other  sources  or 
attenuate  too  quickly. 

The  outgoing  signals  used  here  are  chosen  because  they 
provided  the  best  means  of  studying  the  beacon.  However,  it 
is  likely  these  signals  are  not  the  best  for  operation  in  an 
ocean  or  other  reverberant  environment. 
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Abstract — In  this  paper,  a  method  for  performing  underwater 
acoustic  navigation  using  a  spiral  wave-front  beacon  is  examined. 
A  transducer  designed  to  emit  a  signal  whose  phase  changes  by 
360°  in  one  revolution  can  be  used  in  conjunction  with  a  reference 
signal  to  determine  the  aspect  of  a  remote  receiver  relative  to  the 
beacon.  Experiments  are  conducted  comparing  spiral  wave-front 
beacon  navigation  to  Global  Positioning  System  (GPS)  onboard  an 
unmanned  surface  vehicle.  The  advantages  and  disadvantages  of 
several  outgoing  signals  and  processing  techniques  are  compared. 
The  most  successful  technique  involves  the  use  of  a  phased  array 
projector  utilizing  a  broadband  signal.  Aspect  is  determined  by 
using  a  weighted  mean  over  frequencies.  Sources  of  error  for  each 
of  the  techniques  are  also  examined. 

Index  Terms — Acoustic  devices,  acoustic  navigation,  phased 
array,  underwater  navigation. 


I.  Introduction 

UNDERWATER  navigation  and  localization  techniques  are 
limited  by  attenuation  of  electromagnetic  waves,  espe¬ 
cially  in  saline  or  turbid  water.  Thus,  navigation  aboard  sub¬ 
merged  vessels  uses  a  combination  of  magnetic,  inertial,  and 
acoustic  techniques.  Magnetic  and  gyroscopic  compasses  pro¬ 
vide  only  heading  information.  Inertial  techniques  are  subject 
to  drift  and  must  be  periodically  calibrated  with  absolute  posi¬ 
tion  and  heading  information  [1],  [2].  This  calibration  involves 
surfacing  to  obtain  a  reference  position  from  satellite  signals  or 
using  magnetic  and  acoustic  navigation  information  as  a  refer¬ 
ence.  Thus,  acoustic  navigation  remains  an  important  feature 
in  the  overall  navigation  picture.  Current  acoustic  systems  em¬ 
ploy  arrays,  either  as  projectors  for  long-baseline  navigation  [3], 
or  as  receivers  for  short-baseline  navigation  [4],  [5]  and  ultra¬ 
short-baseline  navigation  [6].  The  arrival  times  or  phases  of  the 
different  elements  are  then  compared  to  determine  positional 
information  using  simple  geometry  [4],  [7],  [8].  In  addition, 
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acoustic  depth  finding  or  sonar  imaging  may  be  compared  with 
maps  to  coordinate  position  and  heading  [9],  [10].  Recent  re¬ 
search  has  focused  on  the  fusion  of  these  different  approaches 
and  the  sharing  of  positional  information  among  vehicles  to  give 
more  accurate  results  [11]. 

The  authors  suggest  an  approach  based  on  a  beacon  capable 
of  producing  a  spiral  wave  front  [12],  [13].  The  beacon  con¬ 
sists  of  two  cylindrical  transducers,  lying  in  the  x-y  plane,  that 
radiate  sound  from  their  outer  edges.  One  transducer,  a  “refer¬ 
ence  source,”  produces  a  wave  front  with  a  cross  section  in  the 
x-y  plane  that  forms  concentric  circles.  A  second  transducer,  a 
“spiral  source,”  produces  a  wave  front  with  a  cross  section  that 
forms  a  linear  spiral.  Consider  a  stable  platform,  whose  position 
is  known,  on  which  the  spiral  beacon  is  fixed  at  a  known  aspect. 
The  beacon  then  sends  a  signal  to  a  remote  receiver  mounted 
to  a  moving  platform,  such  as  an  autonomous  underwater  ve¬ 
hicle  (AUV).  The  AUV  processes  the  signal  to  determine  the 
phase  difference  between  the  sources  that  corresponds  to  the 
AUV’s  aspect  relative  to  that  of  the  beacon.  Range  can  be  deter¬ 
mined  using  time-of-flight  techniques  [14],  [15]  and  depth  from 
a  pressure  sensor.  This  information  can  be  combined  to  give  the 
vehicle’s  position  in  the  3-D  ocean  environment.  In  addition, 
the  navigational  information  gleaned  from  the  spiral  beacon  can 
be  fused  with  navigational  information  from  other  sources  such 
as  Doppler  velocimeters,  gyroscopes,  and  compasses  to  give  a 
more  complete  and  accurate  localization  than  any  of  the  datum 
in  isolation.  This  example  scenario  provides  a  reference  point 
for  the  reader  when  considering  the  experimental  results.  How¬ 
ever,  it  should  be  noted  that  this  is  merely  an  example  con¬ 
figuration.  Various  other  configurations  are  possible,  including 
mounting  the  beacon  onto  a  lead  vehicle,  distributing  a  network 
of  beacons  on  the  seafloor,  or  short-range  vehicle  docking  using 
high-frequency  signals. 

Generating  a  phase  difference  with  the  beacon,  rather  than 
with  an  array  on  the  vehicle  as  with  baseline  techniques,  has 
immediate  advantages.  Most  importantly,  the  remote  platform 
requires  only  a  single  hydrophone  to  determine  its  orientation. 
Often  AUVs  are  equipped  with  hydrophones  for  communica¬ 
tion  or  sonar  systems,  and  these  may  be  utilized  to  receive  the 
beacon  signal.  Also,  any  number  of  vehicles  have  access  to  the 
same  navigational  signal.  In  addition,  it  is  shown  here  that  the 
signal  processing  is  straightforward  and  is,  thus,  easily  adapted 
to  an  onboard  system. 

We  will  also  show  that  the  spiral  wave-front  approach  is  ro¬ 
bust  in  multipath  environments  where  reflected  or  channeled 
rays  are  present.  The  relative  phase  difference  between  the  ref¬ 
erence  and  spiral  sources  is  preserved.  This  sets  the  approach 
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Fig.  1.  (a)  Physical  spiral  wave-front  beacon  showing  both  the  circular  ref¬ 
erence  element  and  the  spiral  element.  When  deployed,  the  spiral  element  is 
located  below  the  reference,  (b)  Phased  spiral  wave-front  beacon.  Below  the 
polyurethane  coating  are  12  evenly  spaced  individual  elements. 


apart  from  baseline  techniques  where  multipath  returns  must  be 
separated  from  the  direct  returns. 

The  concept  of  a  spiral  wavefield  has  been  previously  in¬ 
vestigated  by  Ceperley  [16].  Theoretical  groundwork  for  the 
spiral  wave-front  beacon  is  presented  in  an  introductory  paper 
[12],  and  the  reader  is  referred  there  for  a  complete  discussion 
of  the  concept.  Methods  for  generating  the  spiral  wave  front 
include  a  single  element  transducer  in  the  shape  of  a  linear 
spiral,  a  “physically  induced”  spiral  wave  front,  or  an  array  of 
phased  elements,  a  “phased”  spiral  wave  front.  Two  beacons 
were  developed  utilizing  each  of  these  methods,  and  they  will 
be  referred  to  here  as  the  “physical  beacon”  and  the  “phased 
beacon,”  respectively.  Photographs  of  each  of  the  beacons  are 
shown  in  Fig.  1,  and  detailed  descriptions  and  source  levels  are 
given  in  [17].  Due  to  improper  wiring  of  the  elements  for  the 
phased  beacon,  two  elements  are  shorted  together.  Therefore, 
even  though  there  are  in  total  1 6  elements  available,  neighboring 
pairs  are  driven  together  for  consistency,  giving  in  total  eight 
channels.  As  can  be  seen  in  [17,  Sec.  V  and  Fig.  1],  the  phase 
varies  fairly  smoothly  between  elements  due  to  the  superposi¬ 
tion  of  each  element’s  beam  patterns.  The  propagation  of  a  spiral 
wave-front  signal  in  a  generalized  ocean  channel  is  developed 
in  [1 8].  This  paper  examines  a  set  of  experiments  utilizing  spiral 
wave-front  beacons  for  navigation  aboard  an  unmanned  surface 
vehicle  (USV). 


II.  Beacon  Output  Signals  and  Separation  Techniques 

To  observe  a  phase  difference  between  the  reference  signal 
and  the  spiral  signal,  the  signals  must  be  separated  in  either  time 
or  frequency.  Several  techniques  of  doing  so  were  investigated 
and  a  single  set  of  data  for  each  technique  is  presented  here. 
Three  techniques  utilize  the  physically  induced  spiral  beacon 
and  two  utilize  the  phased  beacon.  The  outgoing  signals  consist 
of  a  series  of  one  or  more  bursts.  In  the  case  of  the  physical 
beacon,  the  spiral  and  reference  transducers  can  send  these 
bursts  out  simultaneously  or  one  at  a  time.  For  the  phased 
beacon  techniques,  we  send  out  first  a  reference  burst  followed 
by  a  spiral  burst,  where  the  spiral  burst  is  generated  by  phase 
delay.  All  of  the  bursts  are  windowed  with  cosine-squared 
shoulders  extending  10%  of  the  pulse  length.  The  reference 
and  spiral  output  signals  of  each  technique  are  shown  in  Fig.  2. 


Time  (ms) 


Fig.  2.  Normalized  output  voltages  for  each  of  the  five  separation  techniques. 
The  reference  signals  are  shown  as  solid  lines  and  the  spiral  signal  is  shown  as 
dashed  lines,  (a)  Separation  in  time  using  the  physical  beacon,  (b)  Separation 
by  frequency  using  the  physical  beacon,  (c)  Frequency-dependent  interference 
using  the  physical  beacon,  (d)  Separation  in  time  using  the  phased  beacon,  (e) 
Separation  in  time  of  a  broadband  signal  using  the  phased  beacon. 
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Fig.  3 .  Normalized  received  voltages  for  each  of  the  five  separation  techniques 
whose  outgoing  signals  are  shown  in  Fig.  2.  Note  that  each  of  these  represen¬ 
tative  arrivals  contains  a  certain  amount  of  reverberation,  primarily  in  the  form 
of  a  strong  surface  reflection.  In  (a),  (b),  and  (c)  the  reflection  does  not  overlap 
the  direct  path  arrival;  in  (d)  and  (e),  however,  it  does. 


Fig.  3  shows  a  representative  arrival  for  each  of  the  techniques 
at  the  vehicle’s  hydrophone  position. 
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Analysis  is  done  offline  after  the  experiments  are  completed 
using  a  desktop  personal  computer  (PC).  For  all  the  techniques, 
the  signal  is  first  bandpass  filtered  to  isolate  the  frequencies  of 
interest.  It  is  then  cross  correlated  with  a  reference  signal  and  its 
envelope  calculated  using  the  Hilbert  transform.  The  first  max¬ 
imum  of  the  envelope  is  used  to  identify  the  location  of  the  first 
burst,  in  the  case  where  there  is  more  than  one  burst.  This  loca¬ 
tion  is  then  used  to  isolate  the  first  arrival  within  a  window  only 
slightly  larger  than  that  of  the  outgoing  burst.  The  windows  for 
later  arrivals  are  determined  by  using  the  known  time  delay  be¬ 
tween  the  bursts.  In  some  cases,  the  computer  records  a  transient 
sound  that  is  not  a  beacon  transmission  (see  Section  III.).  To  rec¬ 
ognize  these  spurious  triggers,  the  normalized  cross  correlation 
between  the  reference  and  spiral  window  is  calculated.  If  this 
value  is  below  0.85,  the  record  is  rejected.  For  techniques  with 
only  a  single  burst,  the  return  can  be  cross  correlated  against  the 
outgoing  signal.  The  Fourier  transform  of  each  window  is  used 
to  determine  the  phase  or  amplitude  required  to  compute  the  as¬ 
pect. 

Each  technique  is  described  in  the  following  paragraphs, 
along  with  the  method  of  aspect  determination.  For  conve¬ 
nience,  throughout  the  paper,  the  techniques  will  be  referred 
to  by  their  letter  below,  i.e.,  (a)  for  the  technique  described  in 
Section  II-A. 

A.  Separation  in  Time  Using  the  Physical  Beacon 

The  most  straightforward  separation  technique  can  be  seen  in 
Fig.  2(a).  Three  successive  tone  bursts  are  sent  from  the  physical 
beacon,  seen  in  Fig.  1(a).  First,  the  reference  transducer  is  acti¬ 
vated,  then  the  spiral,  then  both  together.  Aspect  can  be  deter¬ 
mined  by  using  the  amplitudes  of  each  of  the  three  bursts  or  by 
using  the  phase  difference  between  the  first  two  directly,  as  de¬ 
scribed  in  [17].  The  results  presented  here  use  the  latter  method 
as  it  does  not  suffer  from  left-right  ambiguity  as  the  amplitude 
method  does. 

B.  Separation  in  Frequency  Using  the  Physical  Beacon 

The  arrivals  can  also  be  separated  in  frequency  as  well  as 
time.  For  this  technique,  the  reference  transducer  is  driven  at 
twice  the  frequency  as  that  of  the  spiral  transducer  at  the  same 
time  and  for  the  same  duration.  Aspect  is  determined  using  the 
phase  of  the  received  signal  and  subtracting  the  phase  at  the 
spiral  frequency  from  the  phase  at  the  reference  frequency. 

There  is  difficulty  with  this  technique.  The  phase  difference 
must  be  independent  of  the  window  placement,  which  can  be 
difficult  to  determine  in  a  reverberant  environment.  Indepen¬ 
dence  can  be  maintained  by  using  simple  whole  number  ratios 
between  the  reference  and  spiral  frequencies,  such  as  2:1,  as  is 
done  here.  Other  ratios  lead  to  a  phase  difference  that  appears  as 
a  function  of  window  position.  An  additional  problem  arises  if 
this  technique  is  run  using  the  phased  beacon;  the  angular  side- 
lobe  pattern  is  different  at  different  frequencies  thus  sidelobe 
deviations  are  not  canceled  out  as  they  are  for  same  frequency 
techniques.  For  these  reasons,  this  technique  was  not  successful 
with  the  phased  beacon  and  those  results  are  not  shown  in  this 
paper. 


Fig.  4.  Comparison  of  arrival  spectra  at  various  angles  using  the  frequency-de- 
pendent  interference  technique.  The  arrival’s  spectra  exhibit  aspect  dependence 
due  to  the  superposition  between  the  reference  and  spiral  signals.  Phase  differ¬ 
ence  is  determined  by  comparison  between  received  spectra  (solid  lines)  and 
synthesized  spectra  (dashed  lines). 

C.  Frequency-Dependent  Interference 
Using  the  Physical  Beacon 

The  physical  beacon  is  capable  of  sending  out  the  same 
broadband  signal  from  both  the  reference  transducer  and  the 
spiral  transducer  simultaneously.  Because  there  is  a  path  length 
difference  due  to  the  physical  spiral  that  is  a  function  of  aspect, 
the  frequency  response  changes  with  aspect  due  to  interference 
between  signals.  To  determine  aspect,  the  frequency  response 
of  the  received  signal  is  compared  to  the  frequency  response  of 
synthesized  signals.  Fig.  4  shows  the  frequency  response  of  the 
arrivals  compared  to  that  of  matching  synthesized  signals  for 
three  different  angles.  This  technique  has  left-right  ambiguity. 

D.  Separation  in  Time  Using  the  Phased  Beacon 

Since  each  element  of  the  phased  beacon  can  be  driven  sepa¬ 
rately,  it  is  possible  to  generate  an  outgoing  signal  whose  phase 
advances  by  2i r  radians  across  all  frequencies.  For  an  iV-ele- 
ment  circular  transducer  with  elements  numbered  0  to  N  —  1, 
start  with  a  reference  burst  signal  ro(t)  as  a  function  of  time. 
Then,  we  adopt  the  complex  analytic  signal  as  defined  by 

r  (t)  =  r0(f)  +  iH(r0(t))  (1) 

where  H  (/(£))  denotes  the  Hilbert  transform  of  the  real-valued 
function  f(t).  Then,  the  output  of  each  element  sn(t)  becomes 

sn(t)  =  Re[r  (t)ei2™'N\.  (2) 

In  this  technique,  ro(£)  is  a  simple  tone  burst  with  cosine 
squared  shoulders  extending  1 0%  of  the  pulse  length. 
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At  the  receive  end,  this  technique  is  identical  to  that  using 
the  technique  in  Section  II-A.  To  prevent  reverberation  from  the 
reference  signal  from  interfering  with  the  spiral  signal,  the  time 
between  the  bursts  is  increased  for  the  experiments  utilizing  the 
phased  beacon.  Also,  there  is  no  mixed  third  burst,  because  the 
amplitude  aspect  determination  method  had  not  proven  as  ef¬ 
fective  in  the  earlier  experiments. 


E.  Separation  in  Time  of  a  Broadband  Burst 
Using  the  Phased  Beacon 

Because  the  method  of  driving  the  elements  described  above 
gives  a  phase  shift  across  all  frequencies,  there  is  no  need  to  use 
a  simple  tone  burst.  Thus,  (1)  and  (2)  can  be  used  to  generate  the 
array  output  functions  for  any  finite  pulse.  With  this  technique, 
a  burst  with  a  linearly  swept  frequency  response  is  used.  The 
phase  shift  at  each  frequency  can  then  be  considered. 

Upon  receipt,  the  aspect  is  reported  as  a  weighted  mean  of  the 
phase  differences  across  all  the  frequencies.  The  received  refer¬ 
ence  and  spiral  bursts  r(t)  and  s(t)  can  be  written  as  Fourier 
transforms  R(u)  and  S(uj).  A  weight  is  assigned  to  each  fre¬ 
quency  as  the  product  of  the  two  amplitudes 

W(lj)  =  \R{lj)\\S(lj)\  (3) 

and  the  complex  product  is  defined  as 

P(u>)  =  R(cu)S*(u;)  (3) 


where  S*(uj)  is  the  complex  conjugate  of  the  Fourier  transform 
of  the  spiral  burst.  Aspect  is  determined  by  first  calculating 


(j)n  =  arccos 


Re[P(»] 

W(v) 


then  removing  the  ambiguity 


0'M 


Im[P(cj)]  <  0 
0,  Im[P(u;)]  =  0 

Im[P(a;)]  >  0. 


(4) 


(5) 


The  final  result  is  a  weighed  mean  of  the  amplitudes  across  all 
frequencies 


LU 


(6) 


where  weights  less  than  10%  of  the  maximum  are  rejected  to 
help  eliminate  the  effect  of  noise.  Other  signal  processing  algo¬ 
rithms  can  be  used  for  aspect  determination  that  protect  against 
phase  jumps  across  frequencies.  However,  this  was  not  neces¬ 
sary  here  as  the  beacon  is  operating  far  from  its  resonance. 

It  will  be  shown  that  this  technique  has  two  main  advantages 
over  single  frequency  techniques.  First,  interference  due  to  mul¬ 
tipath  reverberation,  which  can  cause  signal  degradation  at  a 
single  frequency,  is  minimized.  Second,  phase  deviations  due 
to  sidelobe  overlap  of  the  individual  elements  are  reduced. 

It  will  be  shown  that  results  from  the  technique  presented  in 
this  section  [technique  (e)]  do  not  show  the  same  systematic 


errors  seen  in  results  from  the  other  techniques.  For  this  reason, 
we  will  focus  more  carefully  on  these  results. 


III.  Experiment  and  Procedure 

Experiments  utilizing  the  physically  induced  spiral 
wave-front  beacon  (physical  beacon)  took  place  at  the  Naval 
Surface  Warfare  Center-Panama  City  Division  (NSWC-PCD, 
Panama  City,  FL,  USA)  between  August  10  and  20,  2010, 
at  the  Acoustic  Test  Facility.  The  phased  spiral  wave-front 
beacon  experiments  were  carried  out  at  the  Navy’s  Dodge  Pond 
facility  operated  by  the  Naval  Undersea  Warfare  Center-Divi¬ 
sion  Newport  (NUWC-Newport,  Newport,  RI,  USA)  between 
August  9  and  18,  201 1 .  At  each  location,  runs  are  performed  by 
operating  the  USV  by  remote  control  while  it  receives  acoustic 
signals  from  the  beacon.  Only  one  technique  is  tested  during 
each  run,  techniques  (a)-(c)  at  NSWC-PCD,  and  techniques 
(d)  and  (e)  at  Dodge  Pond. 

In  each  of  the  two  sets  of  experiments,  the  beacon  is  fixed 
to  a  floating  dock  at  a  set  depth  below  the  water  line.  Directly 
above  the  beacon  is  a  Global  Positioning  System  (GPS)  antenna 
used  to  precisely  determine  the  beacon’s  position.  The  beacon 
outputs  a  short  circular  reference  burst  followed  by  the  same 
burst,  but  with  a  spiral  wave  front.  These  bursts  are  recorded 
by  a  USV  equipped  with  a  hydrophone  affixed  to  a  rod  below 
the  water’s  surface.  Directly  above  the  hydrophone  is  a  second 
GPS  to  determine  the  hydrophone’s  position.  The  USV  is  driven 
by  remote  control  away  from  the  dock.  Postprocessing  of  the 
signals  allows  the  comparison  between  the  aspect  as  determined 
by  the  GPS  and  the  aspect  as  determined  using  the  spiral  wave- 
front  beacon. 

Fig.  5  shows  a  photograph  of  the  USV  used  in  both  sets  of 
experiments.  The  total  vehicle  length  is  1 .5  m.  The  onboard  sys¬ 
tems  are  summarized  on  the  right-hand  side  of  Fig.  6. 

The  USV  has  two  pontoons  and  is  powered  by  two  small 
side-by-side  trolling  motors  at  the  rear  of  the  boat.  Each  motor 
is  controlled  separately  by  remote  control  allowing  the  vehicle 
to  be  driven  at  a  maximum  range  of  approximately  200  m. 

A  rod  is  attached  to  the  front  of  the  boat  extending  above 
and  below  the  surface  of  the  water.  A  hydrophone  is  attached 
to  the  rod  3  m  below  the  surface  of  the  water.  The  hydrophone 
used  for  the  majority  of  the  experiments  is  an  ITC-1089D  man¬ 
ufactured  by  Channel  Technologies  Group  (Santa  Barbara,  CA, 
USA),  with  an  open  circuit  receiving  response  of  approximately 
— 215  dB  re  1  V//iPa  [19].  For  the  phased  beacon  experiments 
at  NUWC-Newport,  an  ITC-6050C  preamplified  hydrophone  is 
used.  This  transducer  has  an  open  circuit  receiving  response  of 
— 157  dB  re  1  V//iPa  and  its  frequency  range  extends  down  to  4 
kHz  [19]. 

The  signal  received  by  the  hydrophone  is  then  amplified 
through  a  preamplifier  and  recorded  on  an  Agilent  U2531A 
USB  data  acquisition  module.  When  a  transient  signal  comes 
above  a  set  threshold  level,  it  is  recorded  to  a  PC.  Along  with 
the  signal,  the  current  GPS  position  and  the  various  experi¬ 
mental  parameters  are  recorded  simultaneously.  As  the  signal 
level  changes  (for  example,  with  range)  the  threshold  level 
changes  to  anticipate  the  next  signal.  Furthermore,  the  data 
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Spiral  Beacon  is  3  m  below 
the  water's  surface  directly 
below  the  GPS  antenna 
on  the  dock 


Hydrophone  is  3  m  below 
the  water's  surface  directly 
under  the  GPS  antenna  on  the  USV 


Fig.  5.  Annotated  photograph  of  the  USV  and  dock  at  NSWC-PCD  showing 
the  position  of  the  GPS  antennae,  the  system  setup,  and  hardware. 


Reference 
Antenna 
(on  shore) 


Fig.  6.  Cartoon  of  the  data  acquisition,  differential  GPS,  and  control  systems 
onboard  the  vehicle  and  at  the  beacon  position.  The  acoustic  data  and  GPS  co¬ 
ordinates  are  collected  and  synchronized  at  either  end  by  PCs. 


acquisition  module’s  voltage  range  can  shift  dynamically  to  try 
to  maintain  the  maximum  number  of  bits. 

At  the  top  of  the  rod  holding  the  hydrophone,  above  the  water, 
is  a  Novatel  GPS  antenna.  The  GPS  receiver  is  connected  di¬ 
rectly  to  the  antenna,  and  also  via  900-MHz  Freewave  radio  link 
to  a  stationary  antenna  located  off  the  dock  at  a  stationary  loca¬ 
tion.  The  stationary  antenna  provides  a  more  accurate  differen¬ 
tial  GPS  position. 


IV.  Results 

Reported  errors  are  calculated  by  subtracting  the  aspect  as 
determined  by  the  recorded  GPS  positions.  The  beacon’s  exact 
orientation  relative  to  true  north  at  the  time  of  each  run  is  not 


known.  For  the  experiments  with  the  physical  beacon,  an  at¬ 
tempt  was  made  to  align  the  beacon  in  such  a  way  that  the 
physical  gap  [on  the  bottom  right  of  Fig.  1(a)]  was  directed 
away  from  the  USV’s  area  of  operation  to  reduce  diffraction 
effects  from  the  sharp  edge.  For  experiments  with  the  phased 
beacon,  the  zero  aspect  position  was  pointed  away  from  the 
end  of  the  dock  facilities  at  Dodge  Pond  (see  Fig.  8),  roughly 
northwest.  Because  of  these  uncertainties,  both  beacons  must 
be  calibrated  to  determine  the  zero  aspect  orientation.  This  is 
handled  by  adding,  for  each  run,  an  angular  offset  such  that  the 
mean  for  that  run  is  zero.  This  also  corrects  for  experimental  un¬ 
certainty  in  the  positioning  of  the  beacon.  Because  of  this,  the 
root  mean  squared  (RMS)  error  is  a  better  overall  measure  when 
comparing  techniques. 


A.  Sources  of  Error 

Vehicle  motor  noise,  vibration  of  the  hydrophone  support  rod, 
and  local  sound  sources  all  introduce  error  into  the  navigation 
results  and  are  largely  nondeterministic.  However,  there  are  two 
deterministic  mechanisms  that  contribute  to  the  overall  error, 
phase  error  intrinsic  to  the  beacons  and  window  shifting  due 
to  the  radial  velocity  of  the  vehicle.  Because  these  mechanisms 
are  deterministic,  overall  error  could  be  significantly  reduced 
by  compensating  for  them.  However,  since  we  are  interested  in 
studying  and  comparing  the  different  techniques,  this  is  not  done 
in  the  following  results.  Rather,  the  reported  results  reflect  the 
aspect  determination  without  correction. 

One  source  of  deterministic  error  is  due  to  phase  error  in 
the  outgoing  signal.  Both  the  physical  beacon  and  the  phased 
beacon  do  not  generate  perfectly  spiral  outgoing  phases,  as  de¬ 
scribed  in  [17].  For  the  physical  beacon,  these  are  caused  by  ir¬ 
regularities  in  manufacturing,  and,  in  the  phased  beacon,  these 
are  caused  by  sidelobe  overlap  of  the  individual  elements.  These 
errors  are  manifest  in  error  versus  aspect  plots,  as  the  error  is  de¬ 
pendent  on  where  in  the  spiral  field  of  the  hydrophone  resides. 
Phase  errors  are  dependent  on  the  drive  frequency  and,  thus,  are 
minimized  by  looking  at  broadband  signals  where  these  effects 
are  averaged  out. 

Another  source  of  deterministic  error,  for  techniques  where 
the  reference  and  spiral  signals  are  displaced  in  time,  arises 
when  the  radial  velocity  of  the  vehicle  is  nonzero.  The  widowing 
algorithm  described  earlier  assumes  that  there  is  a  fixed  time  be¬ 
tween  reception  of  the  reference  burst  and  the  spiral  burst;  how¬ 
ever,  this  is  not  the  case  when  the  vehicle  has  a  nonzero  radial 
velocity  component  vr  with  respect  to  the  beacon  position.  This 
is  because  the  delay  or  advance  of  the  spiral  burst  causes  the 
window  to  be  misplaced.  This  leads  to  a  phase  error  as  follows: 

A<p  =  (360°).  (7) 

c  —  vr 


In  this  result,  /&  is  the  beacon  center  frequency,  A t0  is  the  time 
between  the  centers  of  the  bursts  at  transmission,  and  c  is  the 
sound  speed  in  water.  If  vr  is  small  compared  to  the  sound 
speed,  a  linear  relationship  can  be  obtained  between  the  phase 
and  the  radial  velocity 


A  0 

vr 


fb^to 


(360°). 


(8) 


c 
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Aspect  (degrees)  Range  (m)  Radial  Velocity  (m/s) 

Fig.  7.  Aspect  errors  from  each  of  the  separation  techniques,  (a)  Separation  in  time  using  the  physical  beacon,  (b)  Separation  in  frequency  using  the  physical 
beacon,  (c)  Frequency-dependent  interference  using  the  physical  beacon,  (d)  Separation  in  time  using  the  phased  beacon,  (e)  Separation  in  time  of  a  broadband 
burst  using  the  phased  beacon  correlated  with  aspect,  range,  and  USV  radial  velocity.  Trend  lines  in  the  radial  velocity  plots  are  calculated  using  the  least  squares 
method. 


This  error  could  be  corrected  easily  by  incorporating  (7)  or 
(8)  into  the  overall  navigation  algorithm.  In  addition,  this  effect 
could  be  minimized  by  using  lower  frequency  signals  and/or 
placing  the  bursts  closer  together.  Finally,  vr  could  be  calculated 
by  correlating  a  second  later  reference  burst  with  the  first 

vr  =  c(6t/ At  —  1)  (9) 

where  St  is  the  measured  time  delay  and  At  is  the  expected  time 
delay. 

B.  Comparison  of  Errors  for  Each  Technique 

Typical  vehicle  paths  for  a  data  collection  run  include  regions 
of  varying  speed,  aspect  angle,  and  range.  Since  the  variety  of 
different  runs  is  large,  a  representative  run  for  each  of  the  tech¬ 
niques  is  chosen  for  discussion  here.  Table  I  gives  the  experi¬ 
mental  conditions  for  each  representative  run.  Plots  of  the  nav¬ 
igational  aspect  errors  for  each  of  these  runs  are  shown  corre¬ 
lated  with  true  aspect,  range  from  the  beacon,  and  USV  speed 
in  Fig.  7.  The  independent  variables  are  taken  from  the  onboard 
GPS.  The  paths  taken  by  the  vehicles,  the  reverberation,  and  ex¬ 
perimental  conditions  are  not  the  same  for  each,  so  these  results 
are  not  directly  comparable.  However,  these  plots  demonstrate 
the  important  features  of  spiral  beacon  navigation.  Table  II  sum¬ 
marizes  the  overall  errors  for  each  run  along  with  five  different 
regions  of  the  data  in  technique  (e). 

For  techniques  (a)-(d),  there  is  a  correlation  between  aspect 
dependence  and  error  due  to  the  phase  error  in  the  outgoing 


Aspect  (degrees) 

Fig.  8.  Errors  calculated  from  simulated  signals  using  [12,  eq.  (30)].  Dashed 
lines  are  calculated  using  the  tone  burst  of  technique  (d),  and  solid  lines  are 
calculated  using  the  linear  swept  frequency  pulse  of  technique  (e).  Phases  are 
calculated  using  the  same  algorithms  used  to  analyze  the  experimental  data.  The 
darker  lines  show  the  results  for  a  1 6-element  array,  and  the  lighter  lines  show 
the  results  for  the  8-element  array  used  in  the  experiments. 

signal,  which  is  an  intrinsic  property  of  the  beacons  [17].  There 
is  less  of  a  correlation  for  technique  (e),  because  the  sidelobe 
effects  are  minimized  by  averaging  over  frequencies. 

No  clear  dependency  is  seen  with  range  for  any  of  the  tech¬ 
niques;  however,  in  the  NSWC-PCD  facility,  the  range  is  lim¬ 
ited  to  57  m. 
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TABLE  I 

Comparison  of  Techniques 


Technique: 

a. 

b. 

c. 

d. 

e. 

Experiment 

Aug, 

Aug, 

Aug, 

Aug., 

Aug., 

Dates 

2010 

2010 

2010 

2011 

2011 

Beacon  Type 

physical 

physical 

physical 

phased 

phased 

(Fig.  la.) 

(Fig.  la.) 

(Fig.  la.) 

(Fig.  lb.) 

(Fig.  lb.) 

No.  of 

2a 

2a 

2a 

8b 

8b 

Channels 

Tone  Type 

tone 

burst 

tone  burst 

linear 

chirp 

tone 

burst 

linear 

chirp 

Cycles  per 

8 

40,  20 

80 

10 

20 

Burst 

Time 

50.6  ps 

N/A 

N/A 

1.33  ms 

1 .20  ms 

Between 

Bursts 

Separation 

time 

frequency 

none 

time 

time 

Frequency 

79  kHz 

159  kHz, 

79  kHz 

40-158 

kHz 

75  kHz 

25-100 

kHz 

p-p  Drive 

200  V 

200  V 

200  V 

230  V 

230  V 

Voltage 

a  One  channel  for  the  reference  signal  and  one  for  the  spiral. 
b  The  transducer  was  designed  to  have  16  elements,  however  two  of 
them  were  wired  together  (see  [17]).  To  compensate  neighboring  element 
pairs  are  driven  together  for  a  total  of  8  channels. 


TABLE  II 

Comparison  of  Techniques  and  Regions  of  Representative  Runs 


Technique 
and  Region 

No.  of 
Datum 

Max. 

Range 

(m) 

Mean 

Speed 

(m/s) 

Mean 

Error 

(degrees) 

RMS  Error 
(degrees) 

a.  (all  data) 

1018 

59 

0.68 

0 

11.1 

b.  (all  data) 

309 

59 

0.63 

0 

12.8 

c.  (all  data) 

326 

58 

0.61 

0 

11.7 

d.  (all  data) 

739 

136 

0.16 

0 

8.3 

e.  (all  data) 

402a 

123 

0.21 

0 

7.2 

e.  Region  I 

8 

70 

1.3 

28.1 

l.lb 

e.  Region  II 

147 

70 

0.074 

2.8 

2.4 

e.  Region  III 

107 

97 

0.16 

1.5 

3.2 

e.  Region  TV 

89 

123 

0.36 

-5.5 

6.8 

e.  Region  V 

28 

28 

0.33 

-6.7 

5.2 

a  Includes  Data  received  between  Regions 
b  Since  only  8  points  are  considered  and  the  aspect  is  constant 


throughout  this  region  this  RMS  error  value  is  unreliable. 

Errors  due  to  radial  velocity  for  techniques  (b)  and  (c)  are 
not  expected  because  A£0  for  these  techniques  is  0.  Equation 
(8)  predicts  the  trend  lines  for  techniques  (a),  (d),  and  (e)  at  3, 
26,  and  23° /(m/s),  respectively.  Although  there  is  considerable 
noise,  the  trend  lines  for  techniques  (d)  and  (e)  confirm  the  pre¬ 
dicted  relationship.  Disagreements  in  the  trends  for  techniques 
(a)  and  (b)  are  likely  due  to  the  limited  aspect  range  utilized. 
The  negative  relationship  observed  with  technique  (c)  is  due  to 
Doppler  shifting  of  the  frequencies.  An  example  of  radial  ve¬ 
locity  error  in  technique  (e)  will  be  pointed  out  for  discussion 
in  Section  IV-C. 

Overall,  results  from  technique  (e)  display  the  lowest  RMS 
error.  In  Section  IV-C,  this  run  will  be  discussed  in  detail. 

The  experimental  errors  from  these  runs  can  be  compared 
to  the  errors  from  theoretical  work  [12]  and  laboratory  testing 
of  the  transducers  [17].  The  RMS  errors  obtained  from  labora¬ 
tory  testing  were  9.1°  for  the  physical  beacon  and  10.1°  for  the 
phased  beacon  run  with  14  single  elements  and  two  shorted  el¬ 
ements  [17].  The  theoretical  error  calculated  in  [12]  was  based 
on  a  transducer  with  16  elements.  In  addition,  those  calcula¬ 
tions  only  relied  on  the  phase  of  a  single  continuous  frequency. 


-150  -100  -50  0  50  100 


East  (m) 

Fig.  9.  Path  of  the  USV  for  technique  (e)  at  Dodge  Pond  in  August  2011.  The 
dock  facilities  are  shown  with  the  beacon  positioned  at  the  origin  denoted  by 
a  white  circle.  The  USV’s  path  as  determined  by  differential  GPS  is  shown  as 
a  white  line.  The  position  as  calculated  at  each  beacon  arrival  is  shown  as  a 
gray  dot.  Range  is  taken  from  the  GPS  result.  Roman  numerals  denote  regions 
discussed  in  Section  IV-C. 

To  simulate  the  expected  error  for  the  phased  beacon,  a  signal 
can  be  simulated  by  using  [12,  eq.  (30)]  as  a  linear  filter  for 
the  reference  burst  ro(t)  for  the  signals  in  techniques  (d)  and 
(e).  The  analytical  model  used  to  calculate  the  beam  patterns  of 
the  elements  assumes  that  the  array  is  baffled,  but  this  shown  to 
be  a  reasonable  approximation  of  the  unbaffled  case  [12].  The 
reference  burst  can  also  be  simulated  with  the  same  equation 
by  using  N  =  1  elements.  The  aspect  can  then  be  calculated 
by  using  the  same  algorithms  used  to  analyze  the  experimental 
data.  These  errors  are  shown  in  Fig.  8.  As  a  comparison,  the  er¬ 
rors  for  a  16-element  array  are  shown  alongside. 

The  RMS  errors  for  technique  (d)  are  9.5°  and  3.0°  for  the 
8-  and  16-element  transducers,  respectively.  For  technique  (e), 
the  errors  are  9.0°  and  2.6°.  Thus,  theoretically,  the  broadband 
weighting  does  not  contribute  much  to  the  overall  error. 

C.  Results  From  the  Separation  in  Time  of  a  Broadband  Burst 
Experiments  at  Dodge  Pond,  August  2011 

Fig.  9  shows  the  path  of  the  USV  at  the  Dodge  Pond  ex¬ 
periment.  Fig.  10  compares  the  aspect  as  determined  by  both 
the  beacon  arrivals  and  the  GPS  onboard  the  USV  and  the 
beacon  for  comparison.  The  roman  numerals  in  Figs.  9  and  10 
denote  the  different  regions  discussed  in  the  following  narra¬ 
tive.  Table  II  contains  the  errors  in  the  individually  discussed 
regions. 

After  leaving  the  dock,  the  vehicle  moves  directly  away  from 
the  dock  at  its  top  speed  of  1 .29  m/s,  marked  as  Region  I  (arrival 
numbers  5 1-58)  in  Figs.  9  and  10.  In  this  region,  there  is  a  mean 
angular  error  of  28.1°.  This  can  be  seen  as  the  group  of  points 
in  Fig.  7,  technique  (e)  versus  USV  radial  velocity  (the  lower, 
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Fig.  10.  Results  from  technique  (e),  run  at  Dodge  Pond  in  August  2011.  (a) 
Compares  the  aspect  as  calculated  by  the  beacon  arrivals  compared  to  aspect 
from  the  GPS  and  (b)  gives  the  corresponding  ranges.  Roman  numerals  denote 
regions  discussed  in  Section  IV-C. 


Sound  speed  (m/s)  Range  (m) 

Fig.  11.  (a)  Sound-speed  profile  from  the  NUWC-Newport’s  Dodge  Pond  fa¬ 
cility  on  the  day  of  technique  (e)  example  run.  (b)  Calculated  transmission  loss 
for  the  sound-speed  profile  given  in  (a)  for  the  direct  acoustic  path  only  at  75  kHz 
for  a  source  at  3-m  depth,  (c)  Calculated  transmission  loss  for  the  sound-speed 
profile  given  in  (a)  for  the  superposition  of  the  direct  acoustic  path  and  the  sur¬ 
face-reflected  path  at  75  kHz  for  a  source  at  3-m  depth. 


rightmost  plot  in  the  figure)  in  the  upper  right  corner.  The  ex¬ 
pected  error  as  calculated  by  (7)  comes  to  29.5°,  indicating  that 
the  error  in  this  region  is  due  to  the  vehicle’s  radial  velocity  in 
this  region. 

After  the  rapid  transit  in  region  I,  the  vehicle  is  allowed  to 
drift  with  the  wind  for  the  next  150  or  so  arrivals.  In  this  region, 
region  II  (arrival  numbers  72-220),  the  beacon  is  the  most  ac¬ 
curate  as  the  USV  is  moving  slowly,  with  a  mean  speed  of  only 
0.074  m/s,  and  there  is  no  noise  from  the  vehicle  motors. 

In  region  III  (arrival  numbers  228-395),  the  vehicle  is  driven 
slowly  at  a  mean  speed  of  0. 16  m/s.  With  this  increase  in  speed, 
the  RMS  error  also  increases.  The  radial  velocity  in  this  region 


is  sufficiently  low  so  that  this  increase  cannot  be  explained  using 
(7).  It  is  most  likely  due  to  motion  of  the  support  rod  holding  the 
hydrophone  and  vehicle  noise. 

In  region  IV  (arrival  numbers  430-608),  the  vehicle  is  driven 
more  aggressively,  with  a  loop,  dramatic  range  variations,  and 
changes  in  vehicle  speed.  This  leads  to  greater  errors  and  a  larger 
number  of  missed  arrivals  as  vehicle  motion  incorrectly  triggers 
the  data  acquisition  system  sporadically. 

As  the  vehicle  returns  home  in  region  V  (arrivals  621-655), 
the  vehicle  turns  from  a  tangential  motion  to  a  radial  one,  which 
leads  to  a  negative  mean  angular  error  in  this  region. 

D.  Performance  of  the  Beacon  in  a  Real  Acoustic  Channel 

Since  the  reference  and  spiral  signals  travel  through  the 
acoustic  channel  in  roughly  the  same  way  [18],  any  phase 
changes  due  to  propagation  are  experienced  by  both  the 
reference  and  spiral  arrivals.  This  makes  the  spiral  beacon  nav¬ 
igation  method  robust  even  in  the  presence  of  a  sound-speed 
profile  or  reflecting  surfaces.  This  can  be  shown  by  examining 
individual  arrivals  at  various  ranges  in  light  of  the  channel’s 
transmission  loss. 

The  sound-speed  profile  at  Dodge  Pond  as  measured  on  the 
day  of  the  experiment  for  technique  (e)  is  plotted  in  Fig.  11(a). 
Fig.  11(b)  shows  the  transmission  loss  in  the  channel  as  calcu¬ 
lated  using  ray-tracing  techniques  for  the  direct  path  only  for 
a  source  located  at  3-m  depth.  A  shadow  region  beginning  at 
around  60-m  range  extends  down  to  hydrophone  depth,  3  m, 
at  around  80  m.  At  this  range,  there  is  no  direct  path  transmis¬ 
sion  from  the  source.  Fig.  11(c)  shows  the  superposition  of  the 
direct  path  with  the  surface  reflection.  Note  that  there  is  inter¬ 
ference  in  the  signals  throughout  the  field  except  in  the  direct 
path  shadow  region.  Thus,  when  the  hydrophone  is  beyond  80 
m,  or  so,  the  majority  of  the  arriving  signal  is  from  the  surface 
reflection.  Because  the  floor  of  Dodge  Pond  is  muddy  and  pri¬ 
marily  absorptive,  reflections  from  the  floor  are  neglected. 

Fig.  12  shows  the  arrivals  at  various  ranges.  Note  that  the 
times  are  shifted  due  to  the  peak  finding  method  described  in 
Section  II-B.  Since  the  arrivals  are  transient  and  not  continuous 
[as  in  the  transmission  loss  in  Fig.  11(b)],  for  short  ranges  there 
is  no  interference  from  the  reflected  path,  as  it  appears  at  a  much 
later  time.  This  is  seen  in  Fig.  12(a).  At  5.7  m,  only  the  direct 
path  signals  contribute  to  aspect  determination.  The  surface  re¬ 
flection  of  the  reference  burst  does  not  appear  until  after  the  di¬ 
rect  path  spiral  arrival  at  5  ms.  The  arrival  in  Fig.  12(a)  and  (f) 
is  number  27  in  Fig.  10. 

In  Fig.  12(b),  at  28.4-m  range,  the  surface  reflection  for  both 
the  reference  burst  and  the  spiral  burst  are  both  clearly  seen  ar¬ 
riving  about  0.5  ms  after  the  direct  path  return.  In  this  case,  the 
algorithms  discussed  in  Section  II-B  can  easily  pick  out  the  cor¬ 
rect  windows.  The  arrival  in  Fig.  12(b)  is  number  27  in  Fig.  10. 
The  arrival  in  Fig.  12(b)  and  (g)  is  number  621  in  Fig.  10. 

In  Fig.  12(c)  and  (d),  the  direct  path  and  the  reflected  path 
overlap  to  varying  degrees  causing  interference,  which  is  a  func¬ 
tion  of  frequency,  as  seen  in  Fig.  12(h)  and  (i).  Since  the  aspect 
is  averaged  over  all  frequencies  in  this  technique,  as  described 
in  Section  II-B,  the  frequency  dropouts  do  not  result  in  any  no¬ 
ticeable  increase  in  error.  Note  that  since  the  window  finding 
algorithm  finds  the  direct  path,  the  window  cuts  the  burst  short, 
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Fig.  12.  Several  arrivals  of  technique  (e)  at  various  ranges.  Graphs  (a)-(e)  show  the  arrivals  in  time.  The  solid  and  dashed  lines  indicate  the  reference  and  spiral 
windows  used  for  aspect  determination  of  aspect,  respectively.  Graphs  (f)-(j)  are  the  respective  frequency  responses  of  the  reference  (solid)  and  the  spiral  (dashed) 
windows. 


just  considering  the  region  in  which  the  direct  return  lies.  The 
arrival  in  Fig.  12(c)  and  (h)  is  number  200  in  Fig.  10,  and  the 
arrival  in  Fig.  12(d)  and  (i)  is  number  301. 

Finally,  in  Fig.  12(e),  when  the  hydrophone  is  well  into  the 
shadow  zone  of  the  direct  path,  there  is  still  a  strong  response 
from  the  reflected  path.  This  is  sufficient  to  receive  a  good  aspect 
result.  Note  also  that  at  these  longer  ranges  there  is  a  degradation 
of  the  higher  frequencies.  The  arrival  in  Fig.  12(e)  and  (j)  is 
number  535  in  Fig.  10. 

There  are  two  caveats  to  consider.  First,  for  the  case  of  single 
frequency  techniques  such  as  (a),  (b),  and  (d),  the  interference 
between  the  direct  and  reflected  paths  can  cause  severe  signal 
degradation  at  ranges  where  the  surface  reflection  cancels  out 
the  direct  path.  Second,  there  can  be  interference  in  the  direct 
spiral  arrival  by  the  surface-reflected  reference  arrival.  Postpro¬ 
cessing  algorithms  recognize  these  situations  and  reject  the  re¬ 
sults  by  comparing  the  energy  levels  in  the  windows. 

V.  Conclusion 

All  techniques  used  for  spiral  navigation  discussed  here  work 
to  determine  aspect  to  varying  degrees  of  success;  however, 
technique  (e)  gives  the  best  results  and  demonstrates  robustness 
in  a  reverberant  channel.  Its  success  is  due  to  the  use  of  a  broad¬ 
band  signal  where  phase  errors  from  sidelobes  and  reduced 
signal  strength  from  interference  are  averaged  out  over  a  band 
of  frequencies.  While  it  appears  from  Fig.  8  and  the  discussion 
in  Section  IV-B  that  averaging  over  a  broadband  gives  only 
marginal  improvement,  consider  Fig.  12(h)— (j).  A  real  acoustic 
channel  can  make  individual  frequencies  inaccessible  due  to 
multipath  reverberation  and  attenuation.  Using  a  broadband 
signal  helps  prevent  dropouts  and  extends  the  working  area  of 
the  beacon. 


For  all  the  techniques,  there  is  stochastic  error  due  primarily 
to  vehicle  noise.  However,  the  majority  of  error  is  from  de¬ 
terministic  sources  as  indicated  by  the  correlations  shown  in 
Fig.  7.  Two  sources  of  deterministic  error  were  identified  and 
discussed:  first,  an  intrinsic  error  that  goes  with  the  AUV’s  ra¬ 
dial  velocity,  and  second,  deviations  of  the  spiral  wave  front 
from  a  pure  linear  spiral.  In  the  former  case,  corrections  can  be 
made  onboard  the  vehicle  in  the  aspect  determination  software. 
In  the  latter,  better  beacon  transducer  design  and  aspect-depen¬ 
dent  corrections  would  minimize  these  errors. 

It  is  hoped  that  study  of  these  results  and  the  lessons  learned 
will  lead  to  the  development  of  a  fieldable  system  for  under¬ 
water  navigation  on  a  submerged  vehicle. 
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